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A computação simbólica ou álgebra computacional é um ramo da Ciência da 


Computação e da Matemática que trata de cálculos simbólicos que envolvem fa- 
toração de polinômios, resolução de equações algébricas e equações diferenciais, 
operações e cálculos com matrizes, grupos, tensores, entre outros. 

Os cálculos realizados no tratamento simbólico são exatos (isto é, têm precisão 
infinita) em contraste ao correspondente tratamento numérico. Embora, em deter- 
minados problemas, pode-se combinar os dois métodos para gerar fórmulas simbó- 
licas que posteriormente serão tomadas com input em programas numéricos, como 
é, por exemplo, o caso da resolução de equações algébricas. 

Nas últimas décadas muitos sistemas de computação simbólica foram desenvolvi- 
dos. Os aplicativos mais conhecidos são o Axiom, Derive, Macsyma, Maple, Matlab, 
Mathematica, Reduce e Maxima. 

Estes aplicativos são de grande utilidade no processo de ensino e na resolução 
de problemas. São ferramentas que podem ser usadas para tornar o ensino da Ma- 
temática mais experimental, além da abordagem de problemas mais complicados 
com visualização gráfica, o que, consequentemente, permite ao estudante aprender 
experimentando. 

Este livro é um guia especial do aplicativo Maxima. Nossa escolha essencialmente 
considerou o fato do Maxima ser um software livre, de acesso gratuito na Internet, 
análogo ao Mathematica e ao Maple e compatível com o UNIX, LINUX, MAC e Win- 
dows. Por ser um software livre possui as seguintes vantagens: acesso ao código 
fonte, incentivo ao desenvolvimento de um espírito de comunidade, maior segu- 
rança, incentivo ao trabalho colaborativo, criação de patrimônio cultural da huma- 
nidade, entre muitas outras. 

Para mais detalhes consulte o sítio http://maxima.sourceforge.net. 

O objetivo principal deste livro é apresentar, e disponibilzar para um público 
mais amplo, as princinpais ferramentas do aplicativo que a autora usou, e ainda 
usa, em sua prática docente, na Universidade Federal do Pará. Portanto, trata-se de 


uma coletânea de tópicos de Matemática que serão apresentados de forma elemen- 


tar dando-se ênfase ao uso do aplicativo. A maioria destes tópicos foram tratados 
pela autora e seus alunos em disciplinas da graduação e orientação de trabalhos de 
conclusão de curso. 

A autora foi inspirada por vários livros, textos e manuais encontrados livremente 


na internet, os quais, em sua maioria, estão listados na bibliografia do livro. 


O livro foi organizado do seguinte modo: no Capítulo[lJapresentaremos os princi- 
pais comandos básicos do aplicativo Maxima: comandos de manipulação algébricas, 
biblioteca e geração de funções, geração de listas, montagem de gráficos e os prin- 
cipais comandos de programação. No Capítulo [2] usaremos o aplicativo Maxima 
para resolver problemas do Cálculo Diferencial e Integral e no Capítulo[3Jusaremos 
o software na geração e programação de Fractais. No final do texto, anexamos a 


bibliografica consultada. 


Belém, 15 de março de 2016. 


Cristina Vaz 





1.1 Informações preliminares 


Um pouco de história 


Maxima é um dos sistemas CAS mais antigos. Foi criado pelo grupo MAC no MIT, 
na década de 1960, e inicialmente chamava-se Macsyma (project MAC's SYmbolic 
MAnipulator). Macsyma foi desenvolvido inicialmente para os computadores de 
grande escala DEC-PDP-10 que eram usados em várias instituições acadêmicas. 

Na década de 1980 foi desenvolvido para várias plataformas, e uma das novas 
versões foi chamada de Maxima. Em 1982, 0 MIT decidiu comercializar Macsyma e, 
paralelamente, o professor William Schelter da Universidade de Texas continuou a 
desenvolver o Maxima. Na segunda metade da década de 1980 apareceram outros 
sistemas CAS proprietários, por exemplo, Maple e Mathematica. Em 1999, a versão 
proprietária do Macsyma foi vendida e retirada do mercado, incapaz de concorrer 
com os outros sistemas CAS. Em 1998, o professor Schelter obteve autorização do 
DOE (Department of Energy), que tinha os direitos de autor sobre a versão original 
do Macsyma, para distribuir o código fonte do Maxima. 

Após a morte do professor Schelter, em 2001, formou-se um grupo de voluntários 


que continuaram a desenvolver e distribuir o Maxima como software livre. 
Instrunções iniciais 


O Maxima trabalha de modo interativo, isto é, ele exibe uma área de trabalho (cha- 
mada wxMaxima), na qual você deverá digitar os comandos e visualizar as respos- 
tas, como mostra a seguinte figura: 


Arquivo Editar Cell Magma Equações Algebra Cálculo Simplificar Gráfico Numérico 
fyuda 


DM ZISDOlRIO DO 16 








Pronto para entrada do usuário 


Algumas instruções preliminares são necessárias antes de você começar a usar O 


aplicativo. 


* Cada linha de entrada (input) é designada por (% in), com n o número da 
linha de entrada e cada linha de saída (output) é designada por (% on), com n 


o número da linha de saída; 
e Cada linha de comando é finalizada com terminador ponto e vírgula (;) 


* Sempre que terminar de digitar uma linha de comando acione a tecla 


shift+ < enter > ou < enter > para que o Maxima efetue as operações; 


* Todos os comandos do Maxima são escritos com letras minúsculas. O apli- 
cativo faz diferença entre maiúsculas e minúsculas. Por exemplo, sin(x) é 


diferente de Sin(x) ou sin(X); 


* Os comandos do Maxima podem ser usados “clicando-se” diretamente nos res- 


pectivos ícones, no menu da área de trabalho. 


1.2 Operações básicas 


1.2.1 Operadores Aritméticos 


Podemos usar o Maxima como uma calculadora efetuando operações numéricas 


com os seguintes operadores aritméticos: 








A - potência 

— - subtração 

/ E divisão 

* = multiplicação 
+ - adição 




















Exemplo 1.1 Efetue as seguintes operações: 


a)5x8x3; b)3,2+6,32; c)(8+52-6x(1+4) d) 





Solução: 





1.2.2 Resultados aproximados e exatos 


Como o Maxima é um aplicativo simbólico, frequentemente vai gerar resultados 


exatos. Para obtermos resultados aproximados usamos um dos seguintes comandos: 


float(M): - resultado decimal do número M. 


M, numer: - resultado decimal do número M. 


1 5 
Exemplo 1.2 Calcule o valor exato e o valor aproximado de sta 


Solução: 


Maxima 
1/5+5/7; 
32 
% Es 
(%o5 ) 35 
float(1/5+5/7); 


(%06) 0.91428571428571 
1/5+5/7, numer; 


(%07) 0.91428571428571 


O aplicativo Maxima pode efetuar cálculos com precisão muito alta. Por exemplo, 
calcular o valor de x com precisão de 40 dígitos: 


Maxima 
fpprec:40; 

(%08) 40 
bfloat(%pi); 


(%09) 3.141592653589793238462643383279502884197b0 


O modificador bfloat foi usado para obter uma representação no formato “big 
float”. A constante fpprec controla o número de algarismos significativos usados 
no formato “big float”. O valor de 40 para a variável fpprec só tem efeito dentro 
do bloco “bfloat” onde foi usado; fora do bloco, fpprec continua com o seu valor 
habitual de 16. A letra b no fim do resultado representa a ordem de grandeza do 


número; neste caso bO representa um fator de 10º = 1. 


Constantes e funções 


O Maxima possui uma biblioteca de constantes e funções que podem ser usadas 


diretamente. 








Constantes usuais 


pi. - €=3.1415926... 
We - ES271626818,. 
Y%l 1=V-1 




















Destacamos os símbolos inf e minf que significam os símbolos matemáticos +00 


e —oo, respectivamente. 








Funções usuais 
abs(x), sqrt(x) = |, 4% 

exp(x), log(x) - e In(x) 

sin(x), cos(x), tan(x) - sen(x), cos(x), tg(x) 
sec(x), csc(x), cot(x) - sec(x), cossec(x), cotg(x). 




















Destacamos os comandos n! e randon(x) que geram o fatorial de n e um número 


aleatório entre [0,x— 1], respectivamente. 


Exemplo 1.3 a) Calcule o valor aproximado de log(11) e e?2; 


b) Calcule o valor exato de sen(=) 
c) Calcule o fatorial de 30. 





Operadores lógicos e relacionais 


Podemos usar o Maxima para efetuar operações lógicas ou relacionais com os se- 


guintes comandos: 








is(expr) - verifica se a expressão expr é falsa ou verdadeira 
assume(expr) - assume que a expressão expr é verdadeira 
forget(expr) - apaga valor lógico da expressão expr 

and - operador lógico e 

or - operador lógico ou 

=, notequal - igual, diferente 

>= - maior , maior ou igual 

Ce - menor, menor ou igual. 




















Não podemos considerar várias expressões condicionais simultaneamente, como 
por exemplo, 6<7< 38. As desigualdades somente se aplicam aos pares de expres- 
sões. Para a resposta verdadeira, o Maxima usa true e para falsa, false. Quando 


acontece uma indefinição, usará a palavra unknown. 


Exemplo 1.4 Efetuando operações lógicas 





1.3 Definindo funções 


Além das funções da biblioteca do Maxima, podemos também definir (criar) novas 


funções usando o comando atribuidor :=. 




















f(x,y... ,Xn):=expr(x,,x2,..., xn) - define uma função f:R” >R 
Flxr. o xn:=[f(X,...,fm(X)] -— - define uma função f:R" > R” 
kill(f) - remove a definição de f 




















Os nomes f que usamos para definir as funções são símbolos que não devem ser 
confundidos pelo programa com nomes já existentes. Por isso, é bom evitar o mesmo 
nome para diferentes funções. Ao terminar de usar o símbolo f, apague sua defini- 


ção usando o comando kill(f). 
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Exemplo 1.5 Defina as funções dadas abaixo e efetue os cálculos indicados. 


2 
a) fo)= ão 


b) g(x,v,z)=y cos(x)+z 
c) h(t) = (t,t7); 0 valor de h(2). 


; 0 valor de f(2.3); o valor exato de f(1); 


2; o valor exato e aproximado de g(5,3,1); 


Maxima 

Ex) c=sqrt((x 2") (322044): 
x2+1 
3x4+1 





(eo lb) flx)= 
E a: 

(%019) 0.2721057705671 
EO 


1 
(%020) — 
v2 
g(x,y,7):=yxcos(x)+202; 
(%021) g(xy,z):=y+cos(x)+z? 
9(5,3,1); 


(%022) 3x+cos(5)+1 
Att e, toa] 


(%023) h(t):= [t,t?] 
h(2) 


(%024) [2,4] 


Funções definidas por sentenças 


Podemos definir funções dadas por sentenças do tipo 


filx) se x<x<x, 
Flx)=4 lx) se x2<x< xs 


Tal) SE Hg SKE 
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Para isto, vamos usar a estrutura condicional if cond then expra else exprb que 
avaliará a expressão expra se a condição cond for verdadeira; se for falsa avaliará a 


expressão exprb. 
Exemplo 1.6 Defina a função 
x2+1 se x<0 


flx)=4 cos(x) se 0<x<r 


2x+3 se x>% 


e calcule f(0), f(—1), f(7). 


Soluções: 





1.4 Manipulando expressões algébricas 


Podemos fatorar, expandir e simplificar uma expressão algébrica com os seguintes 








comandos: 
expand(expr) - efetua produtos e potências da expressão expr 
partfrac(expr var) - simplifica as frações na variável var de expr 
factor(expr) - fatora a expressão expr 
ratsimp(expr) - simplifica a expressão expr 
radcan(expr) - simplifica expressões expr com radicais. 




















Exemplo 1.7 Simplificando expressões. 


Maxima 
p:(x+5)x(x-3); 


(%027) (x-3)x(x+5) 
GE) 2 
(Xo2B) (xr 1) 


ratsimp(p); 


x2+2x-15 


expand(g); 


(%029) 


(%30) x)+3x2+3x41 
factor(%); 
(OS rr 
expand(p/qg); 
2 
2 
(%032) E E 
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fia) 


DD D>>>>—>—>—— 1 
x2 0 2x-1|5 x 20)x- 15 x 42x 15 


Podemos também simplificar expressões trigonométricas com os seguintes co- 


mandos: 


trigexpand(expr) - 


trigsimp(expr) 


efetua produtos e potências sem simplicar expr 


- usaa relação cos?(x) + sen?(x) = 1 para 


simplificar a expressão trigonométrica expr 


Exemplo 1.8 Simplificando expressões trigonométricas. 


Maxima 
trigsimp(tan(a)); 
sin(a) 


sin(b) 
trigexpand(tan(a+b)); 


(%033) 
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tan(b) + tan(a) 
1-tan(a)tan(b) 
trigsimp(2xsin(x)>2+cos(x)22); 


(%034) 


(%035) sin(x)?+1 


1.5 Resolvendo equações 


Podemos resolver equações algébricas usando o seguinte comando: 


solve(eq=0,x) - resolve a equação eq para a variável x. 


A saída do comando solve é uma lista do tipo [x1,x5,x3,...] com x; (i= 1,2,3...) as 
soluções. Se omitirmos a igualdade, o Maxima assumirá o argumento do comando 
solve como uma equação. 


Exemplo 1.9 Resolva as seguintes equações: 


ajreçoels0 bDr=drids0 cjcsly=n=0, 


Maxima 


solve(x4-xr2-1=0,x); 


Vy5+1 Vy5+1 Vy5 — 1%i 
(%o86) [== == CR O, 
2 V2 V2 
Vy5 — 1%i 
o 
V2 
solve(x"5-4+xx+2=0,x); 
(%037) [0=xº-4xx42] 


solve(cos(x)=x,x); 


(%038) [x=cos(x)] 
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1.5.1 Resolvendo equações numericamente 


Note que, o comando solve não resolve as equações dadas nos itens (b) e (c). Po- 


demos resolvê-las numericamente com os seguintes comandos: 








find root(f(x)=0, x, a, b) - método da bissecção; 
newton(f(x)=0,x, xo, imax) - método de Newton; 
mnewton([f,,..., £,) [x1»..., Xu)[x4,...,x2]) - método de Newton 


para sistema de equações. 




















Para o método da bissecção a função deve mudar de sinal nos extremos do inter- 
valo. Se essa condição não for satisfeita, o método será controlado pelo comando 
find root error. Se o valor deste comando for true ocorreram erros; caso contrário 
será retornado um valor. Para o método de Newton, x, é o chute inicial e o critério 
de parada é abs(expr) < imax. Para usarmos o comando newton devemos carregar o 
pacote newtonl e para usarmos o comando mnewton devemos carregar o comando 


mnewton. 


Exemplo 1.10 Resolva numericamente as equações dos itens (b) e (c) do exemplo 
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Exemplo 1.11 Resolva numericamente o seguinte sistema: 


II 
o 


x, +3 In(x1) — 25 


2x2-x1x9-5x1+1 





1.6 Listas 


1.6.1 Gerando listas 


Uma lista para o Maxima é uma série de objetos entre colchetes e separados por 


vírgulas. As listas são geradas pelos seguintes comandos: 








makelist(expr,i,i,,i) - gera uma lista avaliando expr para 
nEle 
create list(expr,iy, Lj,...,iw L,) - gera uma lista avaliando expr com i; 


associado a cada elemento da lista Lj 
e para cada associação indexa i; à cada 
elemento da lista L5,... 




















Exemplo 1.12 (a) Gere uma lista com os 30 primeiros inteiros. 
(b) Gere uma lista formada pelos elementos x, x?, x?, x? 


(c) Obtenha o 15º elemento da lista do item (a) e o 3º elemento da lista do item 


(b). 
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Soluções: 


Maxima 
Li:makelist(i,i,1,30) 


(4042) | [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25, 
26,27,28,29, 90] 
L2:makelist(x 1, 1,4); 
(%043) [x2,x2,x]] 
L1[15];12[3] 


(%044) 15 
3x? 


Exemplo 1.13 (a) Gere uma lista formada pelos elementos (f(1), g(a)), (f(1),g(b)), 
(F (2), g(a)), (02), g(D)). 

(b) Gere uma lista cujos elementos são os gráficos das funções sen(x), 
sen(x + 1), sen(x + 2), sen(x + 2), sen(x + 4) e desenhe esses gráficos com x E [0,27], 
no mesmo sistema de coordenadas. 


Soluções: 


Maxima 
Crcatestis CO on a ea ba 


(%045)  [[f(1), g(a)) [F(1), g(b)], [F(2), g(a)), [f(2), g(b)]] 
plot2d(makelist(sin(k+x),k,0,4),[x,0,2x%pil); 

















(%046) 
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1.6.2 Manipulando listas 


Para manipularmos listas usaremos os seguintes comandos: 








L [il - i-ésimo elemento da lista L 
apply(opr,L) - aplica o operador opr na lista L 
append(L,, L5,..., L,) - juntaaslista Ly, L5,...,L, 


map(opr, L) - “aplica o operador opr na lista L. 




















As operações aritméticas do Maxima são operadores que podem ser usados em lista 
com o comando apply. Por exemplo, apply(“+”,L) aplica o operador soma na lista 
STE: 

Exemplo 1.14 (a) Calcule 12 +23 +32) 44 +53 4+6º+73 +82 +93 +10? 


(b) Gere uma lista com os elementos a, x2, sen(x) e g(x) e calcule o valor de f nessa 
lista. 


(c) Junte as listas [a,b,c,d,e, f,g,h), [1,2,3], [x2,7,32z,w]. 


Soluções: 
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1.6.3 Vetores & Matrizes 


Observe que listas podem ser consideradas “vetores” e podemos efetuar as opera- 
ções usuais dos vetores com os operadores aritméticos do Maxima. Nesse contexto, o 


produto escalar entre vetores é efetuado com o ponto . 


Exemplo 1.15 Considere os vetores u = (x,y,z)e v = (a,b,c). Efetue u+v, u-—v,2u, 


u.v. 


Soluções: 





Matrizes 


Podemos gerar uma matriz com o Maxima usando o seguinte comando: 








matrix([011,012,- R a [4521,922,. ê aim] s[An1,Ãn2» ; E) & 
gera uma matriz nx m. 




















Podemos usar os operadores aritméticos do Maxima para operar com as matrizes. 


O produto entre matrizes será efetuado com o ponto . 
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Exemplo 1.16 Considere as matrizes 


21 
3 5 1 4 6 1 
AS » B=|1 3| €C= 
-2 02 211 
4 1 


Efetue A+C,A-C,5B,A.B. 


Soluções: 


Maxima 
A:matrix([3,5,1],[-2,0,2]) 
3:05 1 
(%055) A= 
-2 02 
BimatrixQpo tos) 


2 
(2056) B=| 13 


RR] 
Comatrix(|4,6,1],[2,1,1]) 


4 61 
oo) 
251] 
A+C 
Pego. 
(%058) 
Gir 
A-C 
-1 -1 0 
(%o59) 
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1.7 Progressões 


Para usarmos os comandos referentes as progressões carregamos o pacote functs. 








arithmetic(a, d,n) - retornao n-ésimo termo da progressão aritmética 
a,a+d,a+2d,..,a+(n-1)d 


geometric(a, r,n) - retorna o n-ésimo termo da progressão geométrica 
a,ar,ar?,..., art 

geosum(a, r,n) - retorna a soma dos elementos da progressão geométrica 

arithsum(a, d,n)  - retorna a soma dos elementos da progressão aritmética 


de 1 an. Sen for infinito (inf) então a soma será finita 
se e somente se o valor absoluto de r for menor que 1. 




















Exemplo 1.17 (a) Determine o décimo termo da PA (2,8,14,...). 


( 
( 
( 
( 
( 


b) Qual a soma dos dez primeiros termos da PA (2,8,14,...). 

c) Calcule a soma dos n primeiros números ímpares (1,3,5,7,...,2n-1,...). 
d) Calcule o décimo quinto termo da PG (1,3,9,27,...). 

e) Calcule a soma dos dez primeiros termos da PG (1,3,9,27,...). 


f) Calcule a soma de todos os termos da PG (1/2,1/4,1/8,...). 
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Soluções: 


Maxima 
load(functs)$ 
arithmetic(2,6,10) 
(%062) 56 
arithsum(2,6,10) 


(%063) 290 
arithsum(1,2,n) 
(%064) nº? 
geometric(1,3,15) 
(%065) 4782969 
geosum(1,3,10) 
(%066) 29524 
geosum(1/2,1/2,inf) 


(%067) 1 


1.8 Gráficos 


1.8.1 Gráficos 2D 


Para desenharmos o gráfico de uma função f(x) (ou mais funções) ou uma curva 
plana definida(s) no intervalo [a,b], em coordenadas cartesianas, usamos os seguin- 


tes comandos: 


plot2d(f(x),[x,a,b], [opções]) - desenha o gráfico da função 
f:[a,b|>R 





plot2d([discrete,Lx,Ly],[opções)) - desenha retas ligando os 
pontos das listas Lx e Ly 
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plot2d([L,],[x,a,b], [opções]) - desenha o gráfico de uma 


lista L, de funções 


plot2d ([parametric,x(t),y(t), [ta,b]) [opções]) - desenha o traço da curva 
(x(t), v(t)) com t e [a,b] 




















Podemos usar diretamente o comando gráfico plot2d do Maxima acionando a ja- 
nela Gráfico 2D localizada na parte superior da área de trabalho. A opção padrão, 
desenha o gráfico numa janela auxiliar e a opção embutido, desenha o gráfico na 
área de trabalho do Maxima. Também podemos obter o gráfico diretamente na área 
de trabalho do Maxima digitando wx seguida do comando plot2d. 

Todas as opções do comando plot2d são listas, sendo a primeira entrada o nome 
da opção. Por exemplo, [nticks, 20] usa 20 pontos para desenhar o gráfico ou o traço 
da curva. O Maxima representa os desenhos do comando plot2d numa proporção de 
4 para 3 entre os eixos horizontal e vertical, portanto para obter a mesma escala nos 
dois eixos usa-se a opção “set size square". 

O Maxima usa os programas auxiliares Gnuplot, Openmath e XMaxima para de- 
senhar os gráficos. No caso do Gnuplot, usa-se o comando gnuplot preamble para 
implementar as opções do plot2d. 


Exemplo 1.18 Desenhe o gráfico das seguintes funções: 


(a) f(x) = sen(2x) e g(x) = sen(3x) no intervalo [0,27] e no mesmo sistema de coor- 
denadas. 


1 
(b)fix)= 7 NO intervalo [-1, 1], com a variação da imagem de f de —100 a 100. 


(c) Desenhe o círculo de centro (0,0) e raio 1. 

















2o 


Maxima 
plot2d(1/x,[x,-1,1],Iy, -100,100]) 


100 














(%069) 


Maxima 
plot2a(]sqriQlo) -sqrt (do o2 o dos dA by, 
[gnuplot preamble, 'set size ratio 1; 


set zeroaxis;"]); 








(%oT70) 

Exemplo 1.19 Desenhe a curva (2t- r.sen(t),2 — 7 cos(t)) com te [-7/2,7/2]. 
Solução: 

Maxima 


plot2d([parametric,2xt-%pixsin(t),2-%pixcos(t),[t,-X%pi,X%pill, 
[nticks,100]) 
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-%pi*cos(t) 


9 











vhOS nora o 
' 


oo 
o 
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o 
[8 
> 
o 
o 


2Pt%pi'sin(t) 


(%071) 


Exemplo 1.20 Desenhe o diagrama que liga os pontos L, = (3,5,3,3,5,3,3,3,1] e 
L, = (3,3,1,3,3,3,5,3,3) e os pontos [((3,1),(5,3),(3,5),(1,3),(3,1); no mesmo sistema 


cartesiano. 


Solução: 


Maxima 

Ex Ss o o Do o o; 
2072) 55,5, 9/9, 3,9, 41] 

Ly SS spo a aa 


(2079) [8:9,1,3,5,9,9,9,9] 
plot2d([ [discrete,Lx,Ly], [discrete,[[3,1],/5,3], 
[551058], [854]]]]); 


5 
discrete1 
45 discrete2 























(%o74) 


2a 


Exemplo 1.21 Desenhe o diagrama que liga os pontos ((0,0), (4,0), (4, 4)), destacando 


as retas e os pontos. 


Solução: 


discrete 1. — 
discrete? é 





1.8.2 Gráficos 3D 


Para desenharmos o gráfico de uma função f(x,y) ou uma superfície definida no 


retângulo [a,b] x [c,d], em coordenadas cartesianas, usamos os seguintes comandos: 








plot3d(f(x,y),|x,a,b], [y,c,d],[opções]) - desenha o gráfico da função 
f:la,b])x[cd]>R 


contour plot(f(x,y),|x,a,b),|y,c,d],[opções]) - desenha as curvas de níveis de 


f(x,y) comxela,b]Jeyvelcd] 
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plot3d([x(u,v),y(u,v), z(u,v)),x,a,b,y,c,d,[opções]) - 


desenha uma superfície parame trizada por ue v com ue [a,bJevelc,d] 




















1-x2-y? 


Exemplo 1.22 Desenhe o gráfico da função f(x,y) = e no domínio 


(xy)e[-3,3]x[-3,3]. 


Solução: 


%eM-y"2-x12+1) 





Exemplo 1.23 Desenhe as curvas de níveis da função f(x,y) = el-=*-»º no domínio 
(x,y)e[-2,2]x[-2,2]. 


Solução: 
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——SppAfyA 2x2 
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1.8.3 Coordenadas polares 


A localização de um ponto P do plano no sistema de coordenadas polares é deter- 
minada pela distância r de P ao pólo O e um ângulo O orientado, no sentido anti- 
horário, formado entre r e a semi-reta horizontal de origem em O, chamada eixo 
polar. Podemos desenhar gráficos em coordenadas polares com o comando plot2d e 


a opção set polar do seguinte modo: 








plot2d([r(0)],[9,0,,05], [gnuplot preamble,“set polar; opções;”]) - 


desenha o gráfico em coordenadas polares. 




















Atenção: para obter a melhor visualização do gráfico, use a opção [7,y1,V2] para 


ajustar a janela de visualização. 

Exemplo 1.24 Desenhe o gráfico em coordenada polares das seguintes funções: 
a) Círculo: r=1; b) Cardióide: r=1 — cos(0) 
c) Rosácea:r = cos(26); d) Concóide: r=2-—sec(0) 


e) Espiral de Arquimedes: r= 0 
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Solução: 


Maxima 
plot2d([1],[ph,0,2x%pi],[y, -2,2], [gnuplot preamble, 


“sei polar” |); 











ph 
(%o79) 


Observe que o círculo foi desenhado como uma elipse!! Para ajustar a janela de 


visulaização usaremos a opção set size square. 


Maxima 
plot2d([1],[ph,0,2x%pi],[y, -2,2], [gnuplot preamble, 


"set polar;set size square"]); 


15 





1+ 
0.5 
= 0 | 
0.5 
1 











15 — 
15 105 005 115 


ph 


(%080) 
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Nos seguintes gráficos, vamos usar as opções unset border e unset tic remover a 
numeração e a borda a implementação gráfica do Maxima. 
Maxima 


plot2d([1-cos(ph)],[ph,0,2x%pil, ly, aa, 
[gnuplot preamble, "set polar; set size square; 
unset border;unset tic;"]); 


1-cos(t) 
e 
e 


(%081) 


Maxima 


plot2d([cos(2+ph)],[ph,0,2x%pil,[y, -2,2], 
[gnuplot preamble,"set polar; set size square; 
unset border;unset tic;"]); 


(%082) 


Maxima 


plot2d([1-sec(ph)],[ph,0,2x%pil,[y, mena 
[gnuplot preamble, "set polar; set size square; 
unset border;unset tic;"]); 


30 


(%083) 


Maxima 
plot2d( [ph], [ph,0,2x%pil, [y, -8,8], 
[gnuplot preamble, 
"set polar; set size square;unset border; 


unset tic;"]); 





(%084) 


Opções do comando plot 


No que segue, destacamos algumas opções do comando plot2d e plot3d. “Op- 
ções” são listas contendo dois ou três elementos, sendo o primeiro elemento o nome 


da opção e os outros os valores associados a esta opção. 


[xax,bx] - amplitude de x 


[y,ay,by] - amplitude de y 
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[nticks,n] - “nm” pontos usados para montagem do gráfico, tem 
padrão n=10 

[grid,n,n] - malha da montagem de gráfico 3D 

gnuplot preamble - insere comandos de montagem; os comandos devem 


estar entre aspas e separados por ponto e vírgula 


gnuplot pm3d - controla o modo avançado PM3D de montagem 
de gráfico 
unset key - omite a fórmula da função. 




















Exemplo 1.25 Desenhe o gráfico e as curvas de níveis da função z = x? +y? no 


mesmo sistema cartesiano. 


Solução: 





Ja 





1.8.4 Gráficos com draw 


O pacote draw é um módulo externo à biblioteca básica do Maxima que permite 
desenharmos gráficos 2D e 3D com relativa comodidade. Este pacote foi exten- 
samente desenvolvido por Mário Riotorto e para ser executado precisa ser pri- 
meiro carregado na memória do computador com o comando load. O pacote é muito 


versátil e possui vários comandos, destacamos os seguintes: 








draw2d(opções, objetos) - desenha os gráficos 2D 


draw3d(opções, objetos) - desenha os gráficos 3D 




















Existem vários comandos que designam as entradas “objetos gráficos”, destaca- 


mos os seguintes: 


Jo 








Objetos gráficos 2D 





explicit(f(x),x,a,b) - gráfico da função f : [a,b] > R; 
parametric(x(t),y(t),tty,t>) - gráfico da curva (x(t), v(t)) com te [ty,t5] 
implicit(eq,x,x1,X2,/Y 152) - traço da curva (x,y(x)) dada implicitamente 


na equação eg com x E [x1,x2],y E [V1,V2] 


polar(r(0), 0,0,,05) - gráfico de uma função r(0) em coordenadas 
polares, com 0 E [04,05] 


points(Lx,Ly) - pontos das listas Lx e Ly que contém os 
valores de x e y, respectivamente 


region(egq(x,y) < 0,x,a,b;y,c,d,) - desenha regiões definidas por desigualdades 
polygon(Lx,Ly) - -poligno com vértices dados nas listas Lx eLy 


rectangle(p,,p>) - retângulo com vértices oposto dados em 
pi=lpx,py] e pz=[px,py] 


ellipse(x9,70,a,b,01,,05) - elipse com centro em (x9,v0) e divisões 
; E a+b 
no eixo de comprimento ERA 
01, ângulo inicial e O, ângulo final 


vector([x1,71])[x2,V2]) - vetor com extremidades (x1,V1) e (x2,72) 




















Exemplo 1.26 Desenhe os seguintes objetos gráficos: 

(a) f(x) = 2sen?(x)+2 cos(x), x e [=x]; 

(b) (x(t), p(t)) = (cos*(t), 2sen?(t)), te [0,27]; 

(c) -xt+y!=1 comy=vy(x), xe[-4,4]eve[-4,4]; 

(d)r = cos(4t), te [0,27]; 

(e) Desenhe todos os objetos gráficos dados acima no mesmo sistema cartesiano; 


2 


as AME 
(f) A região TR HE <1 com (x,v)e[-1,1]x(0,2]. 
y 
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Soluções: 


Maxima 
load(draw)$ 


Solução: (a) 


Maxima 
gl:explicit(2xsin(x)72+2xcos(x),x,-%pi,%pi); 

(%087) explicit(2 + sin(x)? + 2 +cos(x),x, -%pi, %pi) 
draw2d(nticks=200,color=red,g1) 





(%088) 
Solução: (b) 


Maxima 
g2:parametric(cos(t)23,2xsin(t)23,t,0,2x%pi); 

(%089) parametric(cos(t)?,2 «sin(t)?,t,0,2+ %pi) 
draw2d(nticks=200,color=blue,9g2) 








(%090) 


2 


Solução: (c) 


Maxima 
93: implicit(-x-4+y-4=1,x,-4,4,y,-4,4); 
(%091) implicit(-x! +y! = 1,x,-4,4,y,-4,4) 
draw2d(nticks=200,color=green,93) 





4 

a a 
É o Rá 

0 








(%092) 


Solução: (d) 


Maxima 
g4:polar(cos(4+xt),t,0,2x%pi); 
(%093) polar(cos(4+t)+ 1,t,0,2+%pi) 
draw2d(nticks=200,color=orange,g4) 








(%094) 
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Solução: (e) 


Maxima 
draw2d(nticks=200,color=red,g1,color=blue,g2, 


color=green,9g3,color=orange,9g4) 





a 
E tas 


ã / a 


4 
E E 
2 lo e 
1 

0 

-1 








(%095) 


Solução: (f) 


Maxima 
wxdraw2d( 
title = "Gota" ,nticks=400,fill color=magenta, 
region(4xx"2/y-2+y-2/4<=1,x,-1,1,y,.1,2)); 


Gota 














(%096) 


3% 








explicit(f(x,y),x,a,b,y,c,d) 


parametric(x(t),y(t),z(t),tty,t>) 


implicit(eq,x,x1,X2,) 15) 2,2521522) 


points(Lx,Ly, Lz) 


cylindrical(r(z,0),7,71,22,0,01,05) 


spherical(r(0,),0,0,,05,, PP») 


vector([x1,71,21]/[x2,V2,22]) 


Objetos gráficos 3D 





gráfico da função f : [a,b]x[|cd]>R; 


gráfico da curva (x(t), v(t), z(t)) com 
te [ty to] 


traço da curva (x,y, z(x,y)) dada implici, 
tamente na equação eg com x E [x1,%X2] 
evelviv] 


pontos das listas Lx, Ly e Lz que contém 
os valores de x, y e z, respectivamente 


superfície em coordenadas cilíndricas 
com z€E [21,25]; 0 € [04,05] 


superfície em coordenadas espéricas 
com 0 € [01,05], + e hi, Po] 


vetor com extremidades (x1,71,21) 
e (x2,02,22) 




















parametric surface(x(u,v),y(u,v),z(u,V),u,Uj,U5,V,V1,V2) - gráfico da superfície 
parametrizada (x(u,v),v(u,v),z(u,v)) com u e [u,,u5]; v E [04,05] 











Exemplo 1.27 Desenhe os seguintes objetos gráficos: 


(a) f(x,y) = 2sen(xv), xe [-2,2]; 
(b) p=4cos(P), PE [0,27]; 


(Ce 


2 


comxe[-1,1,vel-1,1Jeze[0,1]; 


(d) Desenhe os objetos gráficos dos itens (a) e (b) no mesmo sistema cartesiano. 
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Soluções: 


Maxima 
load(draw)$ 


Solução: (a) 


Maxima 


grli:explicit(2xsin(xxy),x,-2,2,y,-2,2); 


(%097) explicit(2+sin(x+y),x,—2,2,y,-2,2) 
draw3d( (surface hide=true,color=red,gr1)) 


(4098) 
Solução: (b) 


Maxima 
gr2:spherical(4+cos(phi),theta,0,2x%pi,phi,-X%pi,%pi); 
(%099) spherical(4+cos(phi), theta, 0,2 + %pi, phi, -Y%pi, Y%pi); 
draw3d( (surface hide=true,color=blue,gr2)) 


O =to wo 
Climánsr, 





(%0100) 


So 


Solução: (c) 


Maxima 

gro impiltert(xo2 yo r=z o do liyvo- 1, 152,0,d); 
od amplie 20); 

draw3d( (surface hide=true,color=brown,gr3)) 





(%0102) 


Solução: (d) 


Maxima Solução (d) 


wxdraw3d(surface hide=true,color=red,gr1,color=blue,gr2); 





(%0103) 


Opções dos comandos draw2d e draw3d 
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No que segue, destacamos algumas opções do comando draw2d e draw3d. Em ge- 


ral, “opções” são listas contendo dois ou três elementos, sendo o primeiro elemento 


o nome da opção e os outros os valores associados a esta opção. 








color=nome 


fill color=nome 


grid=true/false 
title=nome 


x(yz)label=nome 





Opções gerais 


x(yz)range=[min,max] - 


cor das linhas 

cor de retângulos, polignos e círculos 
amplitude da coordenada x(yz) 
desenha uma grade,se for “true” 
título principal 


título do eixo x(yz) 























filled func=true/false 


filled func= f 


borde=true/false 


line width=valor 








Opções para objetos 2D 


- se for “true”, preenche a área entre o gráfico da 


função e o eixo x 


- se for “true”, preenche a área entre o gráfico da 


função e o gráfico de f 


- se for “true”, desenha a fronteira de polígonos 


- espessura de linhas 
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Opções para objetos 3D 


surface hide=true/false 


controla a visibilidade de partes ocultas 


enhanced3d=true/false se for “true"desenha colorido 


contour=valor - desenha as curvas de níveis; 


valor = none, base, surface, both, map 


contour levels=n - n curvas de níveis defindas por contour 


contour levels=[x,,h,x5] desenha curvas de níveis, defindas por 


contour, entre x; e x, com passo h 


contour levels=[x,,x2,...) - desenha curvas de níveis ,defindas por 


contour, em x4, X5,... 




















X 


Exemplo 1.28 Desenhe em vermelho a área limitada pelo eixo x e a função e” i para 


xe [0,2]. 


Solução: 
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Exemplo 1.29 (a) Desenhe o gráfico e as curvas de níveis da função z = x? +? no 
mesmo sistema cartesiano; 
(b) Desenhe as curvas de níveis da função z = 3sen(|x|) + |y] 
(c) Desenhe uma pirâmide. 


(d) Desenhe um cilidro de raio 1 e altura 6. 


Soluções: 


Maxima 
load(draw)$ 


Solução: (a) 
Maxima 


draw3d(surface hide=true,contour levels=10,contour=base, 
color=pink,explicit(x-24ys2,x,-2,2,y,-2,2)]; 





(%0105) 
Solução: (b) 
Maxima 


draw3d(surface hide=true,contour=map, 


color=purple,explicit(3xsin(abs(x))+abs(y),x,-4,4,y,-4,4)); 
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(40106) 
Solução: (c) 


Maxima 
opcoes: [view=[96,27],surface hide = true,color = green, 
line width = 3]; 

(%0107) [view = [96,27], surfacepide = true, color = green, line width = 3] 
trianguilo(pi;p2,po) «= mesh(|pi,p2l, |ps,po1); 


(%0108) triangulo(p1,p2,p3):=mesh([p1,p2],[p3,p3]) 
quadrado (p1,p2,p3,p4):-mesh([p1,p2],[p3,p4]) 


(%0109) quadrado(p1,p2,p3,p4):= mesh([p1,p2],[p3, p4]); 
piramide : 
Frrtangule([1,0,0]/ 10,150], 1050/11); 
pritangulo( [4,050], 10510], 10,051], 
triangulo([-1,0,0],[0,1,0],[0,0,1]), 
triangulo([-1,0,0],[0,-1,0],[0,0,1]), 


quadrado( [= 10,0] [0, 1, 01/10,=10],11,050]) 18 
draw3d(opcoes,piramide); 


(%0110) 
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Solução (d) 


cilindro 





1.9 Comandos de programação 


É sempre útil fazermos nossos próprios programas, pois otimiza vários processos 
principalmente processos iterativos. Podemos usar a linguagem de programação do 
aplicativo Maxima para implementar algoritmos e resolver vários problemas. No que 


segue, faremos uma breve introdução aos principais comandos de programação. 








if cond then exprl1 else expr2 - avalia a expressão exprl se a condi- 
ção cond for verdadeira; caso contrá- 


rio avalia expr2 


block ([v,,...,V;m) expr,,... expr,) - avalia expr;,... expr, e retorna o 


valor de expr,, 
return(valor) - comando que antecipa uma saída 


/* comentários */ - comentários 
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var:ig, While cond do (expr(var(i))) -avalia a expressão expr para var(i) 
com i iniciando em iç e finaliza quando 


a condição cond for falsa 


var:ig, unless cond do (expr(var(i))) - avalia a expressão expr para var(i) 
com i iniciando em iç e finaliza quando 


a condição cond for verdadeira 


for var(i) in L1 do expr(var(i)) - avalia a expressão expr para 


var(i) nos elementos da lista L1 




















for var(i): io thru i, step m do (expr(var(i))) - 
avalia a expressão expr para var(i) com i variando de iç à i, com passo m, 


finaliza quando i atinge o valor i, 


for var(i): io while cond step m do (expr(var(i))) - 
avalia a expressão expr para var(i) com i iniciando io, com passo m, 


finaliza quando a condição cond for falsa 


for var(i): io unless cond step m do (expr(var(i))) - 
avalia a expressão expr para var(i) com i iniciando io, com passo m, 


finaliza quando a condição cond for verdadeira 




















As declarações do e for são usadas em processo iterativos. As declarações thru, 
while e unless são critérios de parada para os comandos do e for. O processo ite- 
rativo é executado em ciclos: começa com o valor ig para variável var(i) e avalia a 
expressão expr(var(i)); depois um incremento é adicionado ao parâmetro de itera- 
ção, conforme o passo m (se for o caso). Para m=1, o comando step pode ser omitido. 
Em seguida, avalia novamente a expressão expr(var(i)) até que um critério de pa- 
rada seja satisfeito. Para a declaração for, a saída de finalização do Maxima para um 
processo iterativo é done. Em geral, a “programação” com o Maxima gera uma nova 


função. 
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Exemplo 1.30 Faça um programa para calcular n! e calcule 6! 


Solução: 


Maxima 
fatorial(n):= if n < 1 then 1 else n+fatorial(n-1)$ 
fatorial(6) 


(%0112) 720 
Agora, faremos um programa usando os comandos block e do: 


Maxima 
fatoriali(n):=block ([var],var:1,while n> 1 do (var:nxvar, 


n:n-1),var)$ 
fatorial1(6) 


(%0113) 720 


Exemplo 1.31 (a) Some os dez primeiros números naturais; 


(b) Faça um programa para obter o polinômio 4x4 + 7x2 +9x2+ 10x. 


Solução (a): 


Maxima 

block([s],s:0, 

for i: 1 while i <= 10 do (s: s+i),s); 
(%0114) 55 


Solução (b): 


Maxima 

block([p],p:0, 

for i: 1 thru 4 do 

for j: 1 thru i do (p:p+rixxoj),p); 
(xo ls dx 7 OR OK 
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Exemplo 1.32 Faça um programa para calcular os polinônios de Legendre po(x) = 1, 
pi) =, npn(x)=(2n-1)xpn-(x)— (1-1) pa-a(x). Calcule ps(x). 


Solução: 





1.9.1 Aplicações no Cálculo Numérico 


Nesta seção apresentaremos implementações de certos tópicos do Cálculo Numé- 


rico. 


Exemplo 1.33 (aproximação do x) Faça um programa para calcular o valor de 7, 


com precisão dada, usando a série 


a O | 
n=4(1-5+2->+5"""] 
>» dor 4 





Exemplo 1.34 (norma) Faça um programa para calcular a norma euclidiana do ve- 


tor X= (x1,x2,:--,Xxy), dada por 
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Maxima 
norma(lista):=block([n:length(lista),norma], 


soma:sum(lista[il]”2,i,1,n), norma:float(sgrt(soma)))$ 
Ce l,S,5,1,9,14, 18,15, 17 ,19,21,28,25 6 


norma (x); 


(%0118) 50.64582904840241 


Exemplo 1.35 (Método da Bissecção) Faça um programa para calcular uma raiz da 


equação f(x) = x— cos(x), usando o método da Bissecção. 


Solução: 


Maxima /* Método da Bissecção */ 
bissecao(f,a,b,tol, imax):=(x:a, 
for i:1 thru imax do 
(solx[0]:x, x:float((atb)/2),solx[i]:x, 
erro:float(abs((solx[i]-solx[0])/solx[i])), 
if f(a)*f(solx[i]l) < O then b:solx[i] else a: solx[i], 
if erro <= tol then 
(print(“erro relativo=",erro), 
print("solução com a precisão desejada=",x), kill(f),return(fim)) 
else 
if i >= imax then (print("iteração=",i),print("erro relativo=", 
erro), 
print("solução aproximada:",x), 
kill(f), return(fim)) )) $ 


Maxima /* Função */ 
f(x):= x-cos(x)$ 
*/ Verificação se a função muda de sinal nos extremos do intervalo */ 
Fl IE) 


(%0119) —3.721596543649575 
*/ Raíz aproximada */ 
bissecao(f,-1,2,0.001,20); 
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erro relativo=9.904258831297458410-4 
solução com a precisão desejada=0.739501953125 


(%0120) fim 


Exemplo 1.36 (Método de Newton) Faça um programa para calcular uma raiz da 


equação f(x) = x — cos(x), usando o método de Newton. 


Solução: 


Maxima /* Método de Newton */ 
newtonçã x, Lob imp = (ly df aixl] de diet it) +»), 
for i: 1 thru imax do 
(solx]0]:x, y: evldi), x: float(ix = E(o)/y), solx[i]:x, 
erro:abs((solx[i]-solx[0])/solx[il), 
if erro < tol then (print("erro relativo=",erro), 
print("solução com a precisão desejada:",x), 
kill(f),return(fim))))$ 


Maxima /* Função */ 
f(x):= x-cos(x)$ 
*/ Raíz aproximada */ 
newton(£,0.5,0.0001,10); 
erro relativo=7.6489468503479832 10º 
solução com a precisão desejada=0.739501953125 


(%0121) fim 


Exemplo 1.37 (Eliminação Gaussiana) Faça um programa para resolver o sistema 


abaixo, usando o método de Eliminação gaussiana. 


x + 2x) + 3x3 + 4x4 = 10 
2X Fo ds + 2% + Dm —- 7 
3% + 2ão + sy + 2ã = 6 
4x4 + 3x2 + 2x3 + x = 5 


Solução: 








2.1 Somatório e limite 


Podemos efetuar o cálculo de somatórios e limites com os seguintes comandos: 








sum (f(i), i, 1,n) Ê 5 f(i) 


sum(sum(f(i),i, 1, n)j/l,m) - a F(i,9) 
limit(f(x),x,a,plus) - Jim f(x) 
limit(f(x),x,a,minus) “Jim f(x) 
limit(f(x),x,a,inf) “Jim f(x) 


limit(f(x),x,a,minf) -- lim f(x). 


Xx——o0 




















Destacamos os seguintes comandos de saída und e ind para limite indefinido e 


indefinido, mas limitado. 


Exemplo 2.1 Calcule os seguintes limites: 
sen(x). 
x 2, 


a) lim —L; b)lim 
n5+o n+1 x=50 


c) im cos(x). 
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e 


(%0125) 1 


limit(cos(x),x, inf); 


(%0126) ind 


Exemplo 2.2 Calcule os seguintes somatórios: 
10 n 


a) > (+ 1; bb) no c) 3 t+ 
i=1 ist j=1 i=1 
Maxima 


sum(i23+1,1,1,10); 


(%0127) 3035 


sum(i>2,i,1,n) 
n 


(%0128) Rs 


i=1 
sum(sum(i+j,i,1,n),5j,1,m); 
m n 


(40129) DD i+j 


j=1 =] 


Para obtermos o valor simbólico dos somatório (quando houver) com índice 


máximo não númerico, usamos o comando simpsum para simplificar o resultado. 


Maxima 
sum(i”2,i,1,n), simpsum; 
2n2+3n2+n 
(%0130) —————— 
sum(sum(i+j,i,1,n),j,1,m), simpsum; 


mn? mn+mn mn 


%131 
CEO eua na 5 
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2.2 Derivação e integração 


Para calcularmos a derivada e a integral de uma função f usamos os seguintes 
































comandos: 
diff(f(x),x) of) 
diff(f(x) x,n) of) 
; o” [Of 
diff(f(x,y),x,n,y;m) - xi | |) 
b 
integrate(f(x),x,a,b) - J fíx)dx 
b pd 
integrate(integrate(f(x),x,a,b),y,c,d) - J J f(x, v)dxdy 




















Exemplo 2.3 (Cálculo direto de derivadas) Calcule as seguintes derivadas: 
o2f 
ovox 


a)f(x)=vx+1; f(x); Db) f(x,y)=sen(xy?); 
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Podemos definir a função derivada usando o comando define ou a dupla aspas, 
como mostra o seguinte exemplo: 


Exemplo 2.4 (Função derivada) Calcule a derivada da função f(x) = cos(x?) + e 


e o valor desta derivada em x = 1. 


Solução: 


Maxima 
f(x) :=cos(x”2)+exp(-x-2)$ 
define(g(x),diff(f(x),x)); 
(%0135) g(x):=-2xsen(x?) — 2x%e 
diFE(E(x),x)S 
dá (xj=" (0) 
df(1) 
(%0136) -2sin(1)-2%e”! 
9(1) 
(%0137) -—2sin(1)-2%e”| 


Podemos calcular o valor da derivada sem precisarmos definir a função derivada: 


Maxima 
Saad en = 


(%0138) -2sin(1)-2%e”! 
ev(%,x=1); 


(%0139) -—2sin(1)-2%e”1 


Operadores apóstrofo e duplo apóstrofo 


Vamos destacar o uso dos operadores apóstrofo " e duplo apóstrofo ”. O ope- 
rador apóstrofo aplicado a um símbolo não efetua as operações, o operador duplo 


apóstrofo efetua as operações. Veja os seguintes exemplos: 


Maxima 
A:200$ 
A-2 
(%0140) 40000 
"As? 
(XolA) AS 
L “ANP 
(%142) 40000 
“(sqrt(A) + 2xA) 
(%0143) 2A+NVA 
Ros 
(%144) 52%2+400 
SER Ea re O qairi(dO a); 


(%0145) =, (cos(x?) + Ke) =-2xsen(x?)- 2x %e 
x 


Exemplo 2.5 (Cálculo direto de integrais) Calcule as seguintes integrais: 


T/3 x 
a) j sen*(x) cos(x)dx; bb) | mê c) j sds. 
0 x In(x)? o 


Solução: 


Maxima 
integrate(cos(x)xsin(x)”4,x,0, %pi/3); 
5/2 


360 
integrate(1/(x+log(x)22),x); 


(%0146) 


(%0147) 





“Jog6s) 
integrate(s23,5,0,x); 
4 
X 
%o 148 — 
(%0148) A 
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2.2.1 Cálculo de uma variável com o Maxima 


Exemplo 2.6 (Reta tangente e reta normal) Calcule a reta tangente e a reta normal 
à função f(x) = x? no ponto (1,1). Desenhe a função, as retas e o ponto no mesmo 


sistema cartesiano. 


Solução: Para calcularmos a reta tangente vamos derivar a função e obter o coefi- 
ciente angular “m” da tangente; o da normal sabendo que é dado por m, = —1/m. 


Portanto, com o Maxima tem-se: 


Maxima 
Mititae ex 
(%0149) 2x 
difx =" 
(40150) dif(x):=2x 
nd (o: 
(%0151) 2 
retal: 1+mx(x-1); 
(%0152) 2(x-1)+1 
retatg(x):=""%; 
(40158) retaterr=2(K Dl 
reta2:1-(1/m)x(x-1); 


x-—1 
Xo 154 1- 
(%0154) 3 


retanormal(x):=""%; 





(%0155) retanormal(x):= 1 —(x— 1)/2 
Desenhando no mesmo sistema cartesiano: 


Maxima 
load(draw)$ 
grfun:explicit(x>2,x,-1.4,2.8); 
(%0156) explicit(x?,x,-1.4,2.8) 
grtg:explicit(retatg(x),x,-1.4,2.8); 
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retas tangente e normal 





Exemplo 2.7 (Integral definida) Usando a definição, calcule as seguintes integrais: 


a) F x? dx; b) E (x =") da; c) F Vxdx. 


1 
Solução: (a) Devemos calcular o limite: 


n 
ô 
lim ) Si(4-%). 
= 


2 


Como a função f(x) = x* é contínua em [0,1], podemos considerar a partição 


i / ê gu cê ; 
uniforme x;-x; | = — e os pontos €; = x; = —. Assim, o limite acima torna-se 
n n 


n . n 
E 1 e 1 E 1 2 
lim -|) —> lim 5 |. 
n>00 nt n no 2 | 4 
i=1 i=1 
n 
Agora, vamos calcular soma dos quadrados dos n primeiros naturais E, 
i=1 





Calculando o limite: 
1 /2nº+3nº+n E (54 1 ' L) 1 
e lelimi-sre- re AlSs: 
n5013 2n 6nº? 3 


Podemos calcular mais diretamente combinando todos os comandos do seguinte 


modo: 





Solução: (b) Devemos calcular o limite: 


n 


lim (E7- EM)(x;—x;1). 
i=1 


Como a função f(x) = x? x? é contínua em [-1,0], podemos considerar a partição 


. 1 1 . Dea : 
uniforme x;-x; | =— e os pontos €, = x;=—1 + —. Assim, o limite acima torna-se 
n n 


N—00 


n . o. . 3 n . 2 . 
lim »|t+5) -(1+5) E- e (1+-) [1-(=1+-)] 
E] , n n n n>00 1 E n n 


Ki 1 5i 4i2 q? 
O un GR O Ed) 
i=1 
di 12 ud Sa da 
= MDA q Ro EE 
1=— 1— 1— = 
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n 
Resta calcularmos soma dos n primeiros naturais ie a soma dos cubos dos n 
i=1 
n 


primeiros naturais E 


i=1 





Portanto, o limite (2.1) torna-se 





= lim [2->(1+5)e5(2+5 
n>00 2 n 3 nn 

lim |2 E ap Ea +25) 
“* nDo 2 3 4 2n 3n 3nº 4n 4n? 
= (2 E = 
o 2 3 4) 12º 


Podemos calcular mais diretamente combinando todos os comandos do seguinte 


modo: 





Solução: (c) Devemos calcular o limite: 


n 

lim VE = Mia) 
i=1 

Como a função f(x) = yx é contínua em [0,1] podemos considerar a partição 


; L 4 o : 
uniforme x;—- x; | =— e os pontos é; = x; = —. Porém, esta escolha não é a mais 
n n 
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adequada, pois não temos uma fórmula para a série 


n 7 
gas 
o nyn 
1=1 
em termos de “n”. Para evitar este problema, vamos escolher a partição não uni- 
E =P + j2 
fas SM E = > ul e os pontos E; = x; = o 
Assim, o limite acima torna-se 


n 2 & n ç & 
. 1 /21-1 , 1/21-1 
jm al TE Jim to ) 
ca ' rd 


2 1 

q: core 2. Eã na 2.4 E 
e a a A 
= 1= 1= 











Sabemos que 


n n 
1=n, 2 2n2+3n?+n 
6 
1=1 i=1 
Logo, 
n 
. o 2/2n2+3n2+n 4a 
ED E 
1= 
= lim (5+5+ : )=5 
“n5o3 n 3n2) 3 


Podemos calcular mais diretamente combinando todos os comandos do seguinte 


modo: 





Exemplo 2.8 Mostre que 





b 2 2 
bº-a 
E z 


Solução: Usando a definição, devemos calcular o limite: 


lim > & pie 
= 


61 


Como a função f(x) = x é contínua em [a,b|] podemos considerar a partição uniforme 
b-a (b-— ai 


X;—X;1= no eos pontos é; =x;=a+ 





. Assim, o limite acima torna-se 


lim : Eilx;— x;-1) = lim : [a meire, 


n=>00 4 — n5o0 / n n 
i=1 j=1 








Calculando com o Maxima obtemos 


Maxima 
limit(sum((a+(b-a)xi/n)x((b-a)/n),i,1,n),n,inf), 
simpsum; 


b? — Eu 
(40166) d 





Outra solução para este exemplo é usar diretamente o comando integrate: 


Maxima 
integrate(x,x,a,b); 


pn a? 


b 
%0167) >—-— 
Lodo iai 


Exemplo 2.9 (Integral indefinida) Desenhe o gráfico da função 


Você reconhece este gráfico? Explique. 


Solução: 


Maxima 
F(x):= integrate(cos(t),t,0,x) 


assume (x>0); 


(%0168) [x>0] 
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08 
0.6 


E | N F 


sin(x) 
=| 











0.2 
04 N / 
-0.6 te / 
-08 e 
-1 
2 3 4 5 6 7 


(%0 169) 


O gráfico é uma senóide, ou seja, é o gráfico da função f(x) = sen(x). 


Exemplo 2.10 (Integral indefinida) Usando o Teorema Fundamental do Cálculo, 


calcule a derivada das seguintes integrais indefinidas: 
3 


a) Flo) = [ VE + id; brt= |, cos(t) dt. 


r/2 


Solução: (a) 


Maxima 
assume (x>0) 
(%0170) [x>0] 
assume (x>0) 


ratsimp(diff(integrate(sgrt(t>2+1),t,0,x),x)); 


(%o171) Vx241 
Solução: (b) 
Maxima 


diff(integrate(cos(t),t,%pi/2,x73),x); 


(%172) 3x2 cos(x?) 
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Exemplo 2.11 (Cálculo de área) Calcule a área limitada pelo gráfico das funções 


f(x)=2-x2 e g(x) = x. Desenhe a região. 


Solução: Primeiro calcularemos os pontos de intersecção entre as funções: 


Maxima 

Ei(x):=xo [2(x]:=2-x0c2b 

[x1,x2]: map('rhs, solve(£fIlx)=12(x))): 
(%0173) [1,-2] 


Agora, desenharemos a região de integração: 


Maxima 
wxdraw2d(title = "região”,fill color = grey, 
filled func = f2,explicit(fl,x,x2,xd), filled fune = False; 
xaxis = true,xtics axis = true,yaxis = true,line width = 2, 
key = "",color=red,explicit(fl,x,-3,3),key = "color = blue, 


explicit(f2,x,-3,3),yrange=[-4,3]); 


região 

















(%0174) 


1 


Como g(x) < f(x) em [-2,1] temos que a área é J f(x)- g(x)dx. Calculando com 
= 
o Maxima obtemos 


Maxima 
integrate(£2-£1,x,-2,1); 


9 
%0175) > 
Con 


64 


Exemplo 2.12 (Cálculo de área) Calcule a área limitada pelo gráfico das funções 


f(x) = sen(x) e g(x) = cos(x). Desenhe a região. 


Solução: O comando solve do aplicativo Max ima não resolve equações dotipo sen(x) = 
5 


eo aaa T r 
cos(x). Para x € [0,27] temos que o seno e o cosseno são iguais para x = A Car 


Deste modo, a região de integração é dada por 


Maxima 
Ellx)=sin(x)5 t2(x) =cos(x)s 
[x1,x2]: [%pi/4,5+%pi/4]; 


me Sá 
%0176 —, — 
(%0176) [5,>5] 
Maxima 


wxdraw2d(title = "região",fill color = grey,filled func = f2, 
explicit(f1,x,x1,x2),filled func = false,xaxis = true, 


xtics axis = true,yaxis = true,line width = 2,key = ; 


color=red, 


explicit(f1,x,0,5),key="",color = blue,explicit(f2,x,0,5), 


yrange=[-2,2] ); 


região 








(%0 177) 


57/4 
Como sen(x) > cos(x) para x E [7/4,57/4] temos que a área é j (sen(x) — 
r/4 
cos(x)) dx. Calculando com o Maxima tem-se 
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Maxima 
integrate(f1-f2,x,%pi/4,5xX%pi/4); 
(eonne) 22 


Exemplo 2.13 (Cálculo de área) Calcule a área limitada pelo gráfico das funções 


f(x)=3x2-x2-10xe g(x) = —-x? + 2x. Desenhe a região. 


Solução: Primeiro calcularemos os pontos de intersecção entre as funções: 


Maxima 

E1(x):=3:x03-x2-10+x8 E2(x):=-x"2+2+xb 

De loxe,xo0 1: map(' hs, solve(tlbo)=f2(a))); 
(%0179) [-2,2,0] 


Agora, desenharemos a região de integração: 


Maxima 
wxdraw2d(title = "região”,fill color = grey, 
filled func = f2,explicit(f1,x,x1,x2), 
filled func = false,xaxis = true,xtics axis = true, 
yaxis = true,line width = 2,key = "f1",color=red, 
explicit(fl;,x,-3,5), Key="£2";color = blue; 
explicit(f2,x,-3,3),yrange=[-10,7] ); 


região 


» Do 


o 





Sos b to 
1 


' 
— 





(%0180) 


Como g(x) < f(x) em [-2,0] e f(x) < g(x) em [0,2] temos que a área é 


0 2. 
J Ff(x)- (x) d+ | g(x) - f(x)dx. 
=: 0 


Calculando com o Maxima temos 





Observação 2.1 Observe que no exemplo a integração de —2 a 2 daria um re- 
sultado incorreto. De fato, 





Exemplo 2.14 (Cálculo de área) Sejam a > 0 e b > 0. Mostre que a área limitada 


pela elipse 
eu 
o 


2 
1 


érab. 


Solução: Calculando com o Maxima obtemos: 
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Exemplo 2.15 (Cálculo de área) Calcule a área limitada pela da curva abaixo e de- 
senhe a região. 
pesa l=a | 


Solução: Primeiro calcularemos os pontos de intersecção entre as funções 


vy=xVl-x2ey,=-xV1-x? 


Maxima 


assume (-1<x and x<1); 


(%0184) [x>-1,x<1] 
flixesqrt(1-xc2)$f82:-xxsart(1-xr2)$ 
[x1,x2,x9]: map( rhs, solve(f |=[2)): 

(%0185) [-1,1,0] 


Agora, desenharemos a região de integração: 


Maxima 
wxdraw2d(title = "região”,fill color = grey, 
filled func = f2,explicit(f1,x,-1,1), 
filled func = false,xaxis = true,xtics axis = true, 
yaxis = true,line width = 2,key = "f1",color=red, 
explicit(fl,x,-1,1), key="E2",color = blue, 
explicit(f2,x,-1,1),yrange=[-1,1]); 


região 








(%0186) 


Como y, <y em[|-1,0])ey; —y> < g(x) em [0,1] temos que a área é 


0 1 
I w-ndr+ | vi —yv2dx. 
—1 0 


Calculando com o Maxima temos 





Nos seguintes exemplos usaremos o comando integrate para resolver integrais 


indicando qual a técnica de integração usada. 
Exemplo 2.16 (Integração por partes) Calcule as seguintes integrais: 
1 
(a) pe In(x)dx; (b) J arcsen(x)dx  (c) e sen(x) dx. 
0 


Solução: Usando o comando integrate obtemos 





Exemplo 2.17 (Integrais trigonométricas) Calcule as seguintes integrais: 


T/3 3 


cos(x) 
n/6 «/sen(x) 


Solução: Usando o comando integrate obtemos 


(a) ) senta cos(x)tdx; (b) dx (c) | costa! dx. 
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Maxima 
integrate(sin(x)”3xcos(x)”4,x); 
(40191) 5 sat cos(x)> 
integrate(cos(x)23/sgrt(sin(x)),x,%ipi/6,%lpi/3); 
e: 
5923/2 = 523/2 
integrate(cos(x)”4,x),trigreduce; 


(%0192) 


sin(4x) + 8sin(2x) + 12x 


%o 193 
(%01983) 32 


Exemplo 2.18 (Substituição trigonométrica) Calcule as seguintes integrais: 


a | 55 6) [E o | E 


= 3 
dx ——— 
x x2+1)3/2 
Solução: Usando o comando integrate obtemos 





Maxima 
integrate(1/(x-2xsgrt(9-x>2)),x); 


A/ E 
oo sa 


9x 
integrate(sgrt(x”2-3)/x,x,sgrt(3),2),ratsimp; 


V3n — 6 











(%0195) — 
integrate(1/(x-2+1)-(3/2),x); 
x 
(%0196) 
x241 


Exemplo 2.19 (Frações parciais) Calcule as seguintes integrais: 
dx 5x2 +20x+6 8x) + 13x 
S——— (b)|=—s—d SS ——>. 
(a) | a ( ) Dr +x (o | (x2 + 2)2 
Solução: Usando o comando integrate obtemos 


Maxima 


integrate(1/(x"2-D+xx+6),x); 


(%0197) -log(x- 2) +log(x — 3) 
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integrate((D+x72+20xx+6) /(x23+2xx024x),x); 


E 
(%0198) =logtx + 1) + 6log(x) - —— 
integrate((D+x73+13xx) /(x72+2)22,x); 
5log(x? + 2) g 


%199 as 
feonalo 2 2x. +4 


2.2.2 Curvas com o Maxima 


Exemplo 2.20 (Curvas) Seja r(t) = 4cos(t)i + 4sen(t)j + tk. Calcule: 
rá Tr 
a) imo: (55 (9 [Rods dd) | noar 
ÉS dt 0 


Desenhe o traço da curva. 


Soluções: 


Maxima 
r(t):=[4+cos(t),4+xsin(t),t]; 


(%0200) r(t):= [4cos(t), 4sin(t),t] 
limit(r(t),t,%pi); 


(%0201) [-4,0,7] 
GL cRtete), E) 


(%0202) [-4sin(t),4cos(t), 1] 
integrate(r(t),t); 


p 

(%0203) [4sin(t),-4 cos(t), E 
integrate(r(t),t,0,%pi); 

2 

(40204) [0,8,5] 


Desenhando a curva: 


Maxima 
load(draw)$ 
wxdraw3d(nticks=200,parametric(4xcos(t),4xsin(t),t,t,-2,20), 
title="Hélice”"); 
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(40205) 


Exemplo 2.21 Considere a hélice 7(t) = 2cos(t)i + 2sen(t)j + tK. Calcule o vetor T 


. 2, . . L . > 
tangente unitário e o vetor normal unitário N a r(t). Desenhe a curva e os vetores. 


> adiada di adB > vt 
Solução: O vetor T' tangente unitário é definido como T = im e vetor normal 
r 
edi o! RR En > T'(t) 
unitário N é definido como T = am: 
WTC 


Calculando o vetor tangente unitário: 


Maxima 
Ejs [2scos (iu), rs int), te]: 

(%0206) r(t):= [2cos(t),2sin(t),t] 
diff(r(t),t); 

(%0207) [-2sin(t),2cos(t),1] 
define(dr(t),%); 

(%0208) dr(t):= [-2sin(t),2cos(t), 1] 
normadr(t):=sqrt(dr(t)[1]-2+dr(t)[2]-2+dr(t)[3]72); 

(%0209) normadr(t):= «/(dr(t))y + (dr(t) + (dr(t))3 
normadr(t); 

(%0210) 4sin(t)2+4 cos(t)2+1 
Tr(t):=trigsimp(dr(t)/normadr(t)); 


o dr(t 
(%0211) TH) := trigsimp( O) 
Tr(t); 


Za 


2 sin(t) 2 cos(t) 


do 








(ope | ] 


al” 


Calculando o vetor normal unitário: 


Maxima 
dE RRCREÇE), LIS 


(40213) pt 


define(dTr(t),%); 


(%0214)  dTr(t) = [2 SA 520] 


normadir(t):=sgrt(dlr(t)[1]-2+dTr(t)|2]-2+dTr(t)[3]2): 














(%215) normadTr(t) := (dTr(t))Z + (dTr(t)3 + (dTr(t))3 
normadTr(t); 


4 sin(t)? E: 4 cos(t)? 





(%0216) 5 
Nr(t):=trigsimp(dTr(t)/normadTr(t)); 

k e dTr(t) 

(%0217) Nr(9 = trigsimp( a) 


Nr(t); 


(%0218) [-cos(t), -sin(t),0] 


Desenhando a curva e os vetores: 


Maxima 
load(draw)$ 
greurvar:parametrictrit)[U) stt)2la(t)sl,t,-2,10)5 
gel :vector( [O 0/0] [Reto Reto 2 Te(0y STS 
grNr:vector([0,0,0], [Nr(0)[1],Nr(0)[2],Nr(0)[3]])$ 
wxdraw3d(nticks=200,color=green,grcurvar,head length=0.1, 


line width=0.2,color=blue,grTr,color=red,grNr); 
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(%0219) 


Exemplo 2.22 Calcule o comprimento de arco da hélice 


At) = bcos(t)i+ bsen(t)j+ VI — b2tk 
27 
Solução: Sabemos que o comprimento de arco é dado por s = J Ir (ol dt. 
0 


Maxima 
r(t):=[bxcos(t),bxsin(t),sqrt(1-b>2)+t]; 
(%0220) r(t):=[2cos(t),2sin(t), V1 — b2t] 
ea o 
(%0221) [-2sin(t),2cos(t), V1 —b2] 
define(dr(t),%); 
(%0222) dr(t):= [-2sin(t), 2 cos(t), V1 — b2] 
normadr(t):=sqrt(dr(t)[1]-2+dr(t)[2]7-2+dr(t) [3]72); 





(%0223) normadr(t) := (dr(t))2 + (dr(t)2 + (dr(t))3 
normadr(t); 

(%0224) 1 
s:integrate(normadr(t),t,0,2x%pi); 


(%0225) 27 
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2.2.3 Cálculo de várias variáveis com o Maxima 


Exemplo 2.23 (Regra da cadeia) Seja f(x,y)=x?y-—-y? com x=sen(t) e y = e'. 


df 
Calcule ap (O) 


Solução: 


Maxima 
E(x,y):=x02*y -yoZ; 

((Z0226) Rs) e 
[x,y]: [sin(t),exp(t)]; 

(%0227) [sin(t), %e!] 
ditt(f(x,y),t); 

(%0228) %elsen(t)? + 2%e! cos(t)sen(t) — 2 %e?! 
define(df(t),%); 

(%229) df(t) := %e! sen(t)? + 2%e! cos(t)sen(t) — 2%e?! 
af (0) 


(%0230) —2 


d 
Exemplo 2.24 (Derivação implícita) Calcule para y(x) dada implicitamente em 


v'+y?-5y-x2+4=0. 


Solução: 


Maxima 
EV SdVC sy ed 
(%0231) F:iy+y?-5y-x)+4; 
ExsaLttlr, 4) 
(%0232) —2x 
Eycdittdroy): 


(osso 
dy:-Fx/Fy 
2x 


4008/11) 
O mo 
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Exemplo 2.25 (Derivada direcional) Calcule a derivada direcional da função f(x,y) = 


> 


3 5 35 + 
3x? — dor em (59) na direção do vetor v = Z i+j 


Solução: Sabemos que D, f(x,y) = Vf(x,y).if. Usaremos o produto escalar, repre- 


sentado pelo ponto . 





Exemplo 2.26 (Integral dupla) Usando a definição, calcule as seguinte integrais 
duplas: 


(a) [ xvdxdy; (b) ff -v)dxdy; (c) Ff, xy? dxdy. 


Solução: (a) Devemos calcular o limite: 


im[E Ds Vj(x;— x;-1)(9j— vr9) 


i=1 j=1 


Como a função f(x,y) = xy é contínua no retângulo [0,1] x [0,1], podemos consi- 


derar a seguinte partição: 


! j 
E 
E Na 


Assim, o limite acima torna-se 


CAP mal) 


1=1 i=l j=1 


Agora, vamos calcular soma dos n primeiros naturais: 





Assim, devemos calcular o limite: 








Também podemos calcular diretamente: 





Solução: (b) Devemos calcular o limite: 


lim DD) Gx-yla- x) (p= pj) | 


i=1 j=1 


Como a função f(x,y) = 3x — y é contínua no retângulo [1,2] x [-1,1], podemos 


considerar a seguinte partição: 


2i j 
da y=1+5. 
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Assim, devemos calcular o seguinte limite: 


NOM 2i TON] dd Dt 2 
pe [EN fofere 1 DEN] meat 2) 





Solução: (c) Devemos calcular o limite: 
m[Doê vi) sexos) 
t=1 3 


Como a função f(x,y) = x)y? é contínua no retângulo [-1,1] x [0,1], podemos 


considerar a seguinte partição: 


2 
Xj=—, pps 
Assim, 
n n 3 5 2 
te [5 (5) (E) 
e nlAnhn 
3 


Agora, vamos calcular soma dos quadrados e dos cubos dos n primeiros naturais 
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nº+2nº+nº? 


%0245 
(%0245) A 


Assim, devemos calcular o seguinte limite: 


masi pino 


i= j=1 


1 
F 2/1 [nº+2nº+nº 1 Rs 2nº+3n?2+n 
=lm—|—|———— — nes] — 
n>o n2 |n3 4 nº 6 


Maxima 
define(pA(n),%0202)$ 
define(pB(n),%0203)$ 
s(n):=(2xpB(n)/n25)x(1-n+4xpA(n) /n22)$ 


limit(expand(s(n)),n,inf); 


1 
%0246) — 
foston 


Calculando diretamente com o Maxima tem-se: 


Maxima 
limit(sum(sum((i/n)23x(-1+2xj/n)22x(2/n22),i,1,n), 


j,1,n),n,inf), simpsum; 


il 
(%0247) — 
6 


Exemplo 2.27 (Cálculo de volume) Calcule o volume do sólido W limitado abaixo 


pelo retângulo [0,1] x [0,1] e acima pelo plano x+v+z=2 (veja figura/2.1). 


Solução: O volume de W é dado pela integral dupla 


151 
v= [| 2=p=pdxdy: 
040 


Aplicando o Teorema de Fubini obtemos 


ny 





Sólido 


Figura 2.1 


Calculando diretamente com o Maxima tem-se 
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Exemplo 2.28 (Cálculo de volume) Calcule o volume do sólido limitado pelo para- 


bolóide z=4-x? -2y? e pelo plano xy. 


Solução: Fazendo z = 0 temos que a região no plano xy é a elipse x? +2y? = 4, como 


mostra a figura (2.2). 
Neste caso, -2<x<2e qi(x)=-(4-x2)/2e qo(x) = V(4- x?)/2, ou seja, 


4-x?2 A] 











Figura 2.2: Parabolóide elíptico 


Portanto, volume é dado por 
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Exemplo 2.29 (Integral dupla) Desenhe a região de integração D e calcule a inte- 


[sd (x pdxdy. 


Solução: Para desenharmos vamos usar o pacote draw. 


gral 


Maxima 
load(draw)$ 
ERES SIE go 


(%0254) 1-x 
te: sogrt( lx); 


(%0255) Vl-x? 
1340; 


(%0256) O 
wxdraw2d(title = "região”,fill color = grey, 
Eulled func = 3 explicit(t2s,- 1,0), 
filled func = false,xaxis = true, 
yaxis = true,line width = 2,key = "",color=blue, 
explicit(£3,x,-1,0), 


key = ,Color = blue, 
explicit(f2,x,-1,0),1711] color = grey, 

filled func = f3,explicit(f1,x,0,1), 

filled func = false,xtics=true, 

ytics=1/2,lin midth = 2,key = "",color=blie, 
explicit(£3,x,D, 1), key = "color = blue, 
explicht(fl,x,0, |), yrange=[=.3, [9 | grid = troe) 


região 














-0.5 0 0.5 


(%0257) 
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Exemplo 2.30 (Invertendo a ordem de integração) Desenhe a região de integração 


lovx 
D e troque a ordem de integração da integral JJ flsv)dy dx. 
0) x3 


Solução: Para desenharmos vamos usar o pacote draw. 


Maxima 
load(draw)$ 
Elexeas 
f2:sgrt(x)$ 
wxdraw2d(title = "região”,fill color = grey, 
filled func = f1l,explicit(f2,x,0,1), 
filled func = false,xaxis = true, 


yaxis = true,line width = 2,key = ,+Color=blue, 
xtics=true,xrange=[-.1,1.1], 
Ytics=1/2 Iiomidith = 2,key = O scolor=bliye:;, 


explicit(f1,x,0,1),key = ,color = blue, 


explicit(f2,x,0,1),yrange=[-.3,1.3],grid = true); 


região 








(%0258) 


1 yv? 
Trocando a ordem de integração: J J fix, v)dxdy. 
0 v1/3 
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Exemplo 2.31 (Invertendo a ordem de integração) Inverta a ordem de integração 
e depois calcule as seguintes integrais: 


1e1 1e1 
(a) JJ. l+vidydx; (b) JJ Vi -x2dxdy 
0) yx y 


1 py? 


Solução: (a) Invertendo a ordem tem-se I+vidxdy. 


0 
Vamos calcular a integral dupla aplicando o Teorema de Fubini: 


Maxima 
integrate(sqrt(1-y>3),x,0,yº2); 
(%0259) y1l-y3 
f(y):=""%; 


(%0260) f(y):=y? 1 -yº 
Usar a substituição: u=1+v?; du =3y?dy 
integrate(f(y),y,0,1); 


» 
%0261) = 
oco 


Desenhando a região de integração: 


Maxima 
EIS 
f2:sgrt(x)$ 
R6:points([0,0],[0,1])$ 
wxdraw(gr2d(title = "região”,fill color = grey, 
Eniladotume = lsesxplicitit 051); 
filled func = false,xaxis = true,yrange=[-.3,1.3],grid = true, 


yaxis = true,line width = 2,key = ,Color=blue, 


xtics=true,xrange=[-.1,1.1], 


ytics=1/2,line width = 2,key = ,color=blue, 


explicit(fl;x,0,l),key = "* 
explicit(f2,x,0,1), 


color=blue,point type=0,points joined=true,R6)); 


,Color = blue, 
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região 








(%0262) 


lx 
Solução: (b) Invertendo a ordem tem-se j j Vi1-x2dydzx. 
o Jo 


Vamos calcular a integral dupla aplicando o Teorema de Fubini: 


Maxima 
integrate(sqgrt(1-x>2),y,0,x); 
(40263) xNVl-x? 
g(x):=""%; 
(%0264) g(x) =xVI-x2 
sar arsumbstituição: xo du Dada. 


integrate(g(5), x, 0/1); 


1 
(%0265) — 
j 
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Desenhando a região de integração: 


Maxima 
£1:0 
Gia: 
$R7:points([1,1],[1,0])$ 
wxdraw(gr2d(title = "região",fill color = grey, 
filled func = f1l,explicit(f2,x,0,1), 
filled func = false,xaxis = true,yrange=[-.3,1.3],grid = true, 


yaxis = true,line width = 2,key = ,Color=blue, 


xtics=true,xrange=[-.1,1.1], 
ytics=1/2,line width = 2,key = "",color=blue, 


explicit(f1,x,0,1),key = ,color = blue, 
explicit(f2,x,0,1), 


color=blue,point type=0,points joined=true,R7)); 


região 








(%0266) 
Exemplo 2.32 (Cálculo de volume) Calcule o volume do sólido no primeiro octante 
limitado pelas superfícies z = 1-y?,x = y?+1 ex = —y?+9. Faça um esboço do sólido. 


Solução: O volume é dado pela integral da função z = 1 —-y? na região D a ser 


determinada. Primeiro, vamos desenhar a região D. 
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Maxima 


(%0267) 


flisgrt(x-1)$ 


f2:sgrt(9-x)$ 
E3:16 


R8:points([2,21],/0,1])$ 


R9:points([8,8],/0,1])$ 
wxdraw2d(title = "região D",fill color = grey, 


filled func = O,explicit(f1,x,1,2), 

Fill color = grey, Eilled func = O explicit(t3,x,2,0), 
filled func = false,xaxis = true, 

yaxis = true,line width = 2,key = "",color=blue, 
explicit(f3,x,2,8), 


key = ,color = blue, 
explicit(f1,x,1,2),fill color = grey, 
filled func = O,explicit(f2,x,8,9), 
filled func = false,xtics=true, 


ytics=1/2,line width = 2,key = 
explicit(f3,x,2,8), 


,Color=green, 


key = ,color = red, explicit(f1,x,1,2), 


key = ,color = blue, 


explicit(f2,x,8,9),yrange=[0,1.2],grid = true); 


região D 
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A região de integração D é composta de 3 regiões. Portanto, a integral será dada 
por 3 integrais: 1, + 15 + 15; 


1r2 198 1 r-v?+9 
n=)) 1-v?dxdy; b=)), 1-vº dxdy; b=)) 1-v?dxdy. 
0 Jy2+ 0 J2 o Jg 


Maxima 
Il:integrate(integrate(1-y>2,x,y22+1,2),y,0,1); 
8 
(%0268) T5 
I2: integrate(integrate(1-y-2,x,2,8),y,0,1); 
(%0269) 4 
I3:integrate(integrate(1-y>2,x,8,-y-2+9),y,0,1); 


8 
%0270) — 
o 


; 76 
Portanto, o volume é V = ——. 


O sólido é a intersecção de dois cilindros, o que não é trivial desenharmos com o 


Maxima. Para desenharmos este sólido faremos uma reparametrização das superfi- 


cies. 

Maxima 
reparametrize(f1,f2,f83,iv,ivO,ivl,dv,dvO,dv1) := 
apply(parametric surface, append(subst([ iv =u, 
dv = (1-v)xsubst([iv=u],dv0) + v x subst([iv=u],dv1)], 
FERIE2 So us O a Or E; 

(%0271) 


wxdraw3d( 

surface hide = true, 

user preamble = "set view equal xyz;”, 
reparametrize(2xcos(t), 2xsin(t), z, t, O, 2x%pi, Zz, 
-2xabs(sin(t)), 2xabs(sin(t))), 

color = red, 

reparametrize(2+cos(t), z, 2xsin(t), t, O, 2x%pi, z, 
-2xabs(sin(t)), 2xabs(sin(t)))) 
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Exemplo 2.33 (Integral em coordenadas polares) Use integral dupla para calcular 


a área limitada pela rosácea r = 4 cos(4t). Desenhe a região de integração. 


Solução: Sabemos que area = IE dxdy = f fi rdr dt. 
si(t 


Seja 4 a área de uma pétala da rosacea (veja Higura [[3) Então, a região de inte- 
gração é dada por 0O<r <2cos(2,t)e -n/4<t<r/4ea área da pétala é 


T/4 2cos(2t) 
- |) dxdy = J rdr dt. 
R —n/4 0 


Figura 2.3: Rosacea 
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Calculado com o Maxima obtemos 


Maxima 
4xintegrate(integrate(r,r,0,2+xcos(2xt)),t,-X%pi/4,%pi/4); 
(Xo273) 2m 


Exemplo 2.34 (Integral em coordenadas polares) Seja R a região anular limitada 
pelos círculos x? +? = 1 ex? +y? = 4. Calcule a integral | (x? +y)dxdy. Desenhe 


R 
a região de integração. 


Solução: Os extremos polares são 1 <r<2e0<t<27,. Faremos a substituição 


x=rcos(t)ey=rsen(t). Assim, integral é dada por: 


Maxima 
dE (EV) Ra EE 
(%o274) f(x, y):= x? SEN 
pesylolrecos (tt) rrsinçã) | 
(%0275) I[rcos(t),rsen(t)] 


integrate(integrate(f(x,y)xr,1r,1,2),t,0,2x%pi); 


15x 
%o276 —— 
(%0276) A 


Desenhando a região de integração: 


Maxima 


anel 














-2-1.5-1.05 0051152 


(%o277) 


Exemplo 2.35 (Integral tripla) Calcule a integral tripla 


F F [ME “(y+22z)dzdydx 


Solução: 





Exemplo 2.36 (Coordenadas cilíndricas) Calcule o volume do elipsóide 
4x2 +47? +22 =16 


Solução: Sabemos que o volume é dado por 


2 (Vá? p2y4-x2-y? 
v-Jav-s) ) | dzdydx. 
S 0 0 0 


Em coordenadas cilíndricas tem-se 


Tr? r2N4-r? 
v=8 | J J rdzdr dt. 
o Jo Jo 





Exemplo 2.37 (Coordenadas cilíndricas) Considere o sólido S limitado pelo para- 


bolóide z = x? + y? e pelo plano z = 4. Calcule 


JJ. + x +92 da dydx 
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Solução: Usando coordendas cilíndricas tem-se 


4 p2n vz 
|), (2 + va? + vêda dpda = | J 1 ridrdtdz. 
S o Jo Jo 


Maxima 
assume (z>=0)$ 
f(x,y,2):=(x"2+y22) esgrt(x"2+y22)$ 
e,y,z]: [e+cos(t),r=sintt),z]b 
integrate(integrate(integrate(trigsimp(f(x,y,z2))xr,r,0, 
SOPEUZ) Ii ED 2Apij 20 4): 
5127 
39 





(%0280) 


Exemplo 2.38 (Coordenadas esféricas) Calcule em coordanadas esféricas a seguinte 


2 (N4-x2 422? 
J J j vzdz dydx 
-2J0 0 


integral: 


Solução: Usando coordenadas esféricas tem-se 


2 pN4-x? pV4-x2-y? 
J j j vzdzdydx = 
-2 90 0 


Tr/2 pm r2 
e j j (psen(p)sen(0)p cos(qp))p? sen(p) dp do dg. 
0 o Jo 


Maxima 
f(x,y,2) :=yx28 
[x,y,Z]: [rho*sen(phi)*cos(theta),rhoxsin(phi)+xsin(theta), 


rhoxcos (phi)]$ 
integrate(integrate(integrate(f(x,y,z)*xrhoc2xsin(phi),rho,0,2), 


theta,0,%pi),phi,0,%pi/2); 


64 
%281) É 
food as 
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Exemplo 2.39 (Coordenadas esféricas) Calcule o volume do sólido S limitado por 


baixo pelo cone z? = x? + y?,z > 0 e por cima pela esfera x? +y2 +27 = 9 


Solução: Usando coordenadas esféricas tem-se 
m/4 p2r 3 
v-), j j p? sen(qp) dp do dq. 
0 o Jo 


integrate(integrate(integrate(rhor2xsin(phi),rho,0,3), 
theta,0,2+x%pi),phi,0,%pi/4); 


Maxima 


Il 
%o02 82 Dol l= == 
aaa =| 


Exemplo 2.40 (Jacobiano em coordenadas polares) Calcule o jacobiano da mudança 


de variável x=rsen(t)ey =r cos(t). 


Solução: 


Maxima 
xlr,t):=recos(t)S 
ver st) =rasin(t)s 
Jpolar:matrix([diff(x(r,t),r),diff(x(r,t),t)],[diff(y(r,t),r), 
diff(y(r,t),t)]); 


cos(t) —rsen(t) 
(%0283) 


sen(t) rcos(t) 
JaconianoPolar:trigsimp(determinant(Jpolar)); 


(%0284) 1 


Exemplo 2.41 (Jacobiano) Seja Ra região limitada porx-2y = 0,x-2y=4,x+v=4 


ex+y=1. Calcule a integral J 3xydxdy. Desenhe a região de integração. 
R 


Solução: Sabemos que 


| sxvaxa, = Do fato notas do: 
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Maxima 
f(x,y)=3*x*y$ 
x(u,v):=(2xu+v)/3$ 
y(u,v):=(u-v)/3$ 
Jama er ne(uares (x Qa vu) direta quo vo vo ares Qu vu), 
diff(y(u,v),v)]); 
1 


3 


% 
(%0285) 1 


Q|HOTN 


3 
JO1:trigsimp(determinant(J1)); 


1 
%0286) — 
Co. 


integrate(integrate(f(x,y)xabs(J01),v,-4,0),u,1,4); 


164 
(%0287) 


Desenhando a região de integração: 


Maxima 
regi:region(x-2xy<0 and x-2+y>-4and x+y<4 and 


x+y>1,x,-1,3,y,0,3)$ 
retaliexplicit(x/2 x,2/3,0/0)4 


reta2:explicit((x+4)/2,x,-2/3,4/3)$ 
reta3:explicit(4-x,x,4/3,8/3)$ 
reta4:explicit(1-x,x,-2/3,2/3)$ 


Maxima 
wxdraw2d(proportional axes=xy,x voxel=40,y voxe1=40,nticks=450, 
line width=2,fill color=dark-salmon, regi,color=brown, 
retal,color=brown,reta?,color=brown, reta3, color=brown, reta4, 


color=blue, point type=/7,point size=1,A1,B1,C1,D1); 
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(%0288) 


Exemplo 2.42 (Campo vetorial conservativo) Esboce o campo vetorial gerado pelo 
z 
7 


gradiente da função potencial f(x,y) = x2y — 5 


Solução: Usaremos o comando ploteq(f(x,y),|x,y)|x,a,b]),[y,c,d]). Para obter o campo 


clique na grade gerada pelo Maxima. 


Maxima 
pioteqir 2 sy ven lvl Say SO 9] | vectors; blues); 








(%0289) 


5a 


Exemplo 2.43 (Integral de linha) Calcule j (x? -y + 3z)ds com C o segmento de 
reta de extremidades (0,0,0) e (1,2,1). 


Solução: Sabemos que a forma parametrizada do segmento de reta é 7 = (t,2t,t) com 


O<t<le ; 
[tismaas= [ remo) 


com f(x,y,z)=x7-v+32z. 


dr 
Em dt 











Maxima 
f(x,y,72):=x02-y+3x26 
Ee zis bt, 2=t,0]5 
dxditt(|x,y,z];t)6 
linha:integrate(f(x,y,z)xsqrt(dx.dx),t,0,1); 


5 
(%0290) 
V6 


Exemplo 2.44 (Trabalho) Calcule o trabalho realizado pelo campo de forças 


> X 1 


sobre uma partícula que se move ao longo de uma hélice dada por r(t) = (cos(t),sen(t), t). 


Solução: 


Maxima 
Elxy,z)=|-xi2, vi2,1M1]4 
Rey 2: [costt) sin(t) ts 
drsdire(fe, voz],t)S 
W:integrate(F(x,y,z).dr,t,0,3x%pi); 


3n 
%0291 — 
(%0291) A 
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Exemplo 2.45 (Rotacional) O campo vetorial F(x, v,2)=(2xy,x2+22,2y2z) é irrota- 


cional? 


Solução: Mostraremos que rot F = 0. Para isso, usaremos o pacote vect. 


Maxima 
load(vect)$ 
F(x,y,2):=[2xx+y,x02+7202,2xyx2]$ 
cueltE (x,y,2))b 
express(%); 

(%0292) (= (+) lex) 


dz dz 


d 20d 
q (272) qolz + eg 
ev(%, diff); 


X 


(%0293) [0,0,0] 


Exemplo 2.46 (Divergente) Determine o divergente em (2,1,-1) do campo vetorial 


F(x,y, EIS Ea ve Zx22,x? y). 
Solução: 
Maxima 

load(vect)$ 


Elxy,2)=|x0D2y0247,;» 24724702; ne2rvIs 
divergente:diff(F(x,y,z)[1],x)+diff(F(x,y,z)[2],y)+ 
drcrinio yo Z) po zm 

(%0294) 3x2y2z 
ev(%,x=2,y=1,72=-1); 


(40295) —12 
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2.3 Desenhando curvas 
Exemplo 2.47 Desenhe as seguintes curvas: 


(a) (x,y)= (4sen(2t),2 cos(2t)); (b) (x,y) = (cos(t) + tsen(t), sen(t) — t cos(t)); 
31 o] 


(c) (xy)=(Jt-1)t+2); (d) 9)-[ Ton 


Solução: (a) 


Maxima 
load(draw)$ 
curval:parametric(4+sin(2xt),2+cos(2xt),t,0,2x%pi); 
(%0296) parametric(4 sin(2t),2 cos(2t),t,0,27) 
draw2d(nticks=100,color=red,curval); 





o o = 
na wo sa an 


' 
- 











' 
ho 


5 
[ou] 
ho 
dade 
o 
—- 
[0] 
co 
s 


(%0297) 


Solução: (b) 


Maxima 
curva2:parametric(cos(t)+t*xsin(t),sin(t)-txcos(t),t,0,4x%pi); 


(%0298) parametric(tsen(t) + cos(t),sin(t) -t cos(t),t,0,4+7) 
draw2d(nticks=100,color=violet,curva2); 
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(%0299) 
Solução: (c) 


Maxima 


curva3:parametric(abs(t-1),t+2,t,0,2); 


(%0300) parametric(Jt— 1|, t+ 2,t,0,2) 


draw2d(nticks=100,color=green,curvaBS); 


4 








ho 





01 02 03 04 05 06 07 08 09 1 


(%0301) 


Solução: (d) 





Maxima 
curva4:parametric(3+xt/(1+t73),3xt22/(1+t23),t,-12,12); 
So 
(%0302) parametric E RT A lo 


draw2d(nticks=250,xrange=[-3,3],yrange=[0,2],color=orange, 
curva4); 


bo 


ho 








o 
ho 
— 
o 
Es 
ho 
co 


(%03083) 


Exemplo 2.48 A ciclóide é a curva descrita por um ponto P sobre um círculo de raio 
“a” rolando ao longo de uma reta num plano. As equações paramétricas da ciclóide 
são 
(x,y) = a(t —- sen(t)),a(l — cos(t)). 
Desenhe a ciclóide (x,y) = 2(t — sen(t)),2(1 — cos(t)) 


Solução: 


Maxima 
load(draw)$ 
cicloide:parametric(2+x(t-sin(t)),2x(1-cos(t)),t,0,8x%pi); 


(%0304) parametric(2 (t — sen(t)),2(1 — cos(t)),t,0,87) 
draw2d(nticks=100,xrange=[0,8+x%pi],yrange=[0,5], 
title="ciclóide”,color=coral,cicloide); 


ciclóide 














0 5 10 15 20 25 


(%0305) 
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Exemplo 2.49 A epiciclóide é a curva descrita por um ponto P sobre um círculo de 
raio “b” que rola do lado de fora de um círculo de raio “a” com a > b. As equações 


paramétricas da epiciclóide são 





x = (a+b)cos(t)-b cos 


Es 


2) 


) 





v = (a+b)sen(t)-bsen 


Ss 


— —"— 
q 
+ 
Ss 


Desenhe a epiciclóide 


1 
É = Hcos(t) -3cos(); 
1 
v = Lsen(t) -3sen(—), 
Solução: 
Maxima 
load(draw)$ 


set draw defaults (user preamble="set size ratio -1",nticks=200, 


dimensions=[300,300],title="epiciclóide",color=red)$ 
epicicloide:parametric(11*xcos(t)-3xcos(11/3xt), 


lixsin(t)-3*sin(11/3+xt),t,0,6x%pi); 


é 11t 11t 
(%0306) parametric(11 cos(t) —3 cos(——), 11 sen(t) — 3sen( >). 0,67) 


wxdraw2d(epicicloide); 


epiciclóide 


(%0307) 
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Exemplo 2.50 A hipociclóide é a curva descrita por um ponto P sobre um círculo 


de raio “b” que rola no interior de um círculo de raio “a” com a > b. As equações 


paramétricas da epiciclóide são 





e (= b)cos(8) + b cos[S 2); 





—b 
v = (= b)sen()- bsen( ) 
Desenhe a hipociclóide 


q 
x = 5cos(t)+3 cos(=); 


v = Ssen(t) -3sen( >) 
Solução: 


Maxima 
load(draw)$ 
set draw defaults(user preamble="set size ratio -1",nticks=200, 


dimensions=[300,300],title="hipociclóide”,color=salmon)$ 
hipocicloide:parametric(bxcos(t)+3x+cos(D/3xt), 


bxsin(t)-3+sin(5/3xt),t,0,6x%pi); 


Su Su 
(%0308) parametric(5 cos(t)+ 3 cos(=), 5sen(t) — 3sen(=)t 0,67) 
wxdraw2d(hipocicloide); 


hipociclóide 


1 1 1 1 
oo EmonmnmEeE oo 
1 


(%0309) 
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Exemplo 2.51 Desenhe as seguintes curvas: 
2 
(a) rosácea: r=cos(5) com 0 e [0,67]; 
(b) Espiral hiperbólica: r = 2/0; 
(c) Estrofóide: r = 2 cos(20) sec(0). 


Solução: (a) 


Maxima 
load(draw)$ 
set draw defaults(user preamble="set size ratio -1",nticks=200, 
dimensions=[350,350],title="rosácea",color=plum)$ 
rosacea: polar(cos(2xth/3),th,0,6x%pi); 


(%0310) polar [cos[2 un 0,67] 


wxdraw2d(rosacea) ; 
rosácea 


0.8 
0.6 
0.4 
0.2 


-0.2 
-04 
-0.6 
-0.6 


-1-0.8-0.60402 0 02040608 1 
(%0311) 
Solução: (b) 


Maxima 


espiral: polar(2/th,th,0,8x%pi); 


E 
(%0312) polar (5th, 0,87) 
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wxdraw2d(nticks=200,title="espiral hiperbólica”,color=turquoise, 
xrange=[-1,3],yrange=[-2,2],espiral); 


espiral hiperbólica 








(%0313) 


Solução: (c) 


Maxima 
estrofoide:polar(2+cos(2xth)xsec(th),th,0,2x%pi); 
(%0314) polar(2 cos(2th) sec(th),th, 0,27) 
wxdraw2d(nticks=200,title="estrofóide”,color=brown, 


xrange=[-2,3],yrange=[-5,5] ,estrofoide); 


estrofóide 








(%0315) 
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2.4 Desenhando superfícies 
Exemplo 2.52 Desenhe o plano x —- 5y+22z = 0. 


Solução: 


Maxima 
plano1:x-D+y+2xz =0; 


(40316) 2z-5y+x=0 
draw3d(enhanced3d = true, implicit(plano1,x,-4,4,y,-4,4, z,-6,6)); 





(%0317) 

x2 v2 2 
Exemplo 2.53 Desenhe o elipsóide Sa ao 1. 
Solução: 
Maxima 


elipsoide1l:x»2/2+y-2/3+772/4=1; 
draw3d(enhanced3d =true, 
implicit(elipsoidel,x,-2,2,y,-2,2,27,-2,2)); 
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Exemplo 2.54 Desenhe a parte do sólido no 1º octante, limitado pelas superfícies 


x=v2+1,x=-y?+9ez=1-y?, 


Solução: Vamos parametrizar as superfícies, considerando O< y<1le0<z<1. 
Assim, 

Sr: x=v+1>v=tx=t+1;0<t<1l;0<z<1-t; 

93º x=-v2+9>5y=t,x=-t"+9;0<t< lO<ez<i=t; 

Sa: el-rei 

Para obtermos a variação da variável “x”, vamos calcular as interseções entre as 
superfícies S; e S; e as superfícies S, e S5. Assim, 

Sins dons SÊ; 

S,NS3:x-z=8>5x=9-t, 

Logo, O<t<1;l+t?<x<9-t2. 

Para desenharmos este sólido usaremos o comando draw3d. Como os extremos 
das superfícies são funções do parâmetro, faremos uma reparematrização para ter- 
mos extremos fixos. 

Primeiro, desenharemos as superfícies do tipo grade e depois do tipo cheia. Para 
melhor vizualização, desenharemos o plano xz e no tipo grade desenharemos cada 


superfície de uma cor diferente. 
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Maxima 
reparametrize(f1,f2,f3,iv,ivO,ivi,dv,dv0O,dv1) := 
apply(parametric surface, append(subst([ iv =u, 
dv = (1-v)xsubst([iv=u],dv0) + v * subst([iv=u],dv1)], 
RESP ss) us vos ivo vo O ld: 

(%0319) out 
wxdraw3d(surface hide = true, 
reparametrize(t”2+1,t,z, t, 0, 1, z,0,1-t72), 
color = red, reparametrize(-t-2+40,t,2, t, 0, 1, z,0,1-t2), 
color = yellow,reparametrize(x,t,1-t>2, t, 0,1, x, 1+t22, 
9=1"2), 


color = green,reparametrize(x,0,7, x, 1,9, 2, 0, 1)) 





(%0320) 
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Maxima 
wxdraw3d(enhanced3d = true, 
reparametrize(t>2+1,t,z, t, 0, 1, z,0,1-t72), 
color = red,reparametrize(-t>2+49,t,z, t, 0, 1, z,0,1-t72), 
color = yellow,reparametrize(x,t,1-t22, t, 0,1, x, 1+t22, 
9-t"2), 


color = green,reparametrize(x,0,7, x, 1,9, 2, 0, 1)) 


S000 
O NEMO 

o9900- 
him o 





(%0321) 


Exemplo 2.55 


(a) Desenhe o cilindro senoidal y = 2 sen(x); 
(b) Desenhe o Hiperbolóide z = x? -y2, (x,y) E [-2,2] x [-2,2] e suas curvas de níveis. 


Solução: (a) 


Maxima 
draw3d(enhanced3d=false,surface hide-=false, 


parametric surface(t,6xsin(t),u,t,0,2x%pi,u,0,2)); 
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(%0322) 


Solução: (b) 


Maxima 
wxdraw3d(enhanced3d=false, 
colorbox=false, surface hide=true, contour=base, 
explicit(xr2-y>2,x,-2,2,y,-2,2)) 





(%0323) 
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Exemplo 2.56 Parabolóide hiperbólico é uma superfície de revolução obtida girando 


a hiperbóle (t,1/t) em torno do eixo x. Uma parametrização desta superfície é 


* 


Ê Ê 


[r cos(v) e) 


coml<t<2e0<v<2r. 
Desenhe o parabolóide usando esta parametrização em torno dos eixos y e z. De- 


senhe as curvas de níveis. 


Solução: 


Maxima 
wxdraw3d(enhanced3d=false,colorbox=false, surface hide=true, 
contour=base,parametric surface 


(u,cos(v)/u,sin(v)/u,u,1,2,v,0,2x%pi)) 











SAS 
CS aa aa 
COZ 4, E aa 
Tiaras 
ar A e a 






(%0324) 


Maxima 
wxdraw3d(enhanced3d=false, 
colorbox=false, surface hide=true, contour=base, 
parametric surface(cos(v)/u,sin(v)/u,u,u,1,2, 
v,0,2x%pi)) 
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Exemplo 2.57 Uma parametrização do Elipsóide é 
(asen(u)cos(v),bsen(u)sen(v),ccos(u)) 
com (u,v) € [0,27]x[0,27]ea,b,ceR. Desenhe o elipsóide e as curvas de níveis. 


Solução: 
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Exemplo 2.58 Uma parametrização do Cone é 
(ucos(v),usen(v),u) 


com (u,v) E [a,b] x [0,2%pi] e a,b E R. Desenhe o cone no formato grade de no 
formado cheio. Desenhe curvas de níveis. 


Solução: 
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(40327) 


Maxima 


SA pIGUICS, 


wxdraw3d (enhanced3d 


=base, 


true, contour 


hide= 


parametric surface(uxcos(v),uxsin(v),u,u,-2,2, 


v,0,2x%pi)) 


false, surface 


colorbox 





(%0328) 
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Exemplo 2.59 O Catenóide é uma superficie de revolução obtida girando-se a cate- 
nária (t,acosh(t/a)) em torno do eixo x, com -b<t<b e a,b > 0. Uma parametriza- 


ção do catenóide é 
(t,(acosh(t/a)) cos(w), (acosh(t/a)) sen(w)) 


com -b<t<be0O<wx<27r. Desenhe o catenóide em torno dos eixos x e z. Desenhe 


curvas de níveis. 


Solução: 


Maxima 
wxdraw3d (enhanced3d=false, 
colorbox=false, surface hide=true, contour=base, 
parametric surface(t,cosh(t)xcos(w),cosh(t)+sin(w), 
t,-1,1,w,0,2x%pi)) 





(%0329) 


Maxima 
wxdraw3d(enhanced3d=true, 
colorbox=false, surface hide=true, contour=base, 
parametric surface(cosh(t)xcos(w),cosh(t)xsin(w), 
t,t,-1,1,w,0,2x%pi)) 
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Exemplo 2.60 O Toro é uma superfície de revolução obtida girando-se um círculo 
de raio r e centro (0,4, 0) em torno do eixo z, com a >r > 0. Uma parametrização do 
Toro é 


((a+rcos(v))cos(u), (a+rcos(v))sen(u), rsen(v). 


Desenhe o Toro nos formatos grade e cheio. No formato grade, desenhe curvas de 


níveis. 


Solução: 
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(%0331) 


Maxima 


(%0332) 


wxdraw3d(enhanced3d=true, 
colorbox=false, surface hide-=true, 
parametric surface((I+cos(v))x*cos(u), 


(l+cos(v))*sin(u),sin(v),u,0,2x%pi,v,0,2x%pi)) 


05 —— 
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Exemplo 2.61 Uma parametrização do Helicóide é 
(ucos(v),usen(v),v) 
coma<u<xbe0O<v<2r. Desenhe o Helicóide nos formatos grade de cheio. 
Desenhe curvas de níveis. 
Solução: 
Maxima 
wxdraw3d (enhanced3d=false, 
colorbox=false, surface hide=true, contour=base, 


parametric surface(uxcos(v),uxsin(v),v,u,-2,2, 
v,0,2x%pi)) 





(%0333) 
Exemplo 2.62 Uma parametrização da Faixa de Móbius é 
((a+usen(v/2))cos(v), (a+usen(v/2))sen(v),ucos(v/2)) 


com -a<u<ae0)O<v<2r,. Desenhe a Faixa de Moóbius. 


Maxima 
wxdraw3d(enhanced3d=false,surface hide=true, 
parametric surface((5+uxsin(v/2))xcos(v), 
(5+uxsin(v/2))xsin(v),uxcos(v/2),u,-1,1,v,0,2x%pi)) 
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2.5 Equações diferenciais ordinárias 


Para resolvermos equações diferenciais ordinárias com o Maxima usamos o se- 


guinte comando: 








ode2(edo, vardep, varindep) 


edo - equação diferencial 
vardep - variável dependente 


varindep - variável independente. 




















Quando o aplicativo resolve a EDO, retorna uma solução (explícita ou implícita) 
para a variável dependente e representa as constantes por “Y%cl", “%c2", etc. Se não 
pode obter a solução, retorna a palavra false seguida, às vezes, de uma mensagem 


de erro. Para obtermos o valor das constantes usamos o seguinte comando: 
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icj(sol, varindep=x9, vardep=vo) 


j=1 - EDO de primeira ordem. 

j=2 - EDO de segunda ordem. 

sol - solução geral da EDO. 

xo - valor inicial da variável independente. 


Vo - valor inicial da variável dependente. 




















Para o caso de problemas de fronteira com EDO de segunda ordem, obtemos o 


valor das constantes usando o seguinte comando: 








bc2(sol, frontx=x;, fronty=y,, frontx=x», fronty=y>) 


sol - solução geral da EDO 
x1, X> - valores de fronteira da variável independente 


V1, V2 - Valores fronteira da variável dependente. 




















Exemplo 2.63 Encontre a solução geral das seguintes equações diferenciais lineares: 
(a) v'+ay=0; 
(b) y'+ap=f(0); 


(o) v'+a(x)y = f(x). 


com a uma constante não nula e f(x) uma função dada. 
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Solução: (a) 


Maxima 


eql:“diff(y,x)+axy=0; 


d 
A =, = (() 
(%0335) Ray 
solt:ode2(eql,y,x); 


[%O396) y— %0c%6€7"* 
Solução: (b) 


Maxima 
depends(f,x); 
(%0337) [f(x)] 
eq2:"diff(y,x)+raxy=f(x); 
d 
(%0338) o = f(x) 
sol2:ode2(eqg2,y,x); 


(%0339) y= toe | oe-t* tg dx + Yc 


Solução: (c) 


Maxima 
depends(f,x); 
(%0340) [f(x] 
depends(a,x); 
(%0341) [a(x)] 
eq03:diff(y,x)+ra(x)xy=f(x); 
(%0342) Cpra(ay= f(0 
so103:0de?2(eq03,y,x); 


) aa - J ata 
(%0343) y=Y%e Jão f(x) dx + %c 
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Exemplo 2.64 Considere a equação 


dy eus 
Gê : 


(a) Ache a solução geral; 


(b) Ache a solução y(x) que satisfaz a condição inicial v(0) = O. Esboce o gráfico. 


Solução: (a) 





Solução: (b) 
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(%0348) 


Exemplo 2.65 Resolva a equação 


dy “ 
PR 


Solução: 
Maxima 
eq04: diff(y,x)2xxxy=x; 


d 
(%0349) O =X 
so104:0de2(eq04,y,x); 


He 
(%0350) p= me [EE sue 


Exemplo 2.66 Considere a seguinte equação 


dy Euaddl 
o AO 


chamada equação de Bernoulli. 


(a) Ache a solução geral; 


(b) Ache a solução que satisfaz v(0) = 1. 
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Solução: (a) 


Maxima 
eq02: diff(y,x)+2xxxy=xxy22; 
d 
(%0351) ay +20 =xy 
sol05:0de?2(eq05,y,x); 


“log(p+2)-log(y) x? 


(%0352) 5 = +%e 


Solução: (b) 


Maxima 
dell sol05, x=0,y=1); 


Jog(y+2)-log(y) x? —log(3) 


(%0353) 
p 2 





Exemplo 2.67 Encontre a solução v(x) do seguinte problema de valor inicial 


(x-Dy?+(y-1xy'=0, 
v(3)=2 


Solução: 


Maxima 
eq: (x-1)xy23+(y-1)xx23x" diff(y,x)=0; 


(%0354) Pp +(e-193=0 
ode? (eg,y,x); 
(%0355) = = he 
peltcos x 24 o) 
27-1 x2+72x-36 
De Mo e 








(%0356) 
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Exemplo 2.68 Encontre a solução geral da equação diferencial de 2a ordem com 
coeficientes constantes dada por 


v"+by'+cy=0. 


Solução: 


Maxima 
eq4: diff(y,x,2)+bx' diff (y,x)+cxy=0; 
(%0357) En bx] +cy=0 
dx? dx 
sol4:0de2(eq4,y,x); 
positive 
No =] 
j 


X 4 = b? 
(%0358) y= %e É | %k1 sin + %k2 cos E) 
sol4:0de2(eq4,y,x); 
negative 


(Nb2-4-b)x (-Vb2-4c-b)x 
(%0359) y=%kl%e 2 + %k2%e 2 


sol4:0de?(eg4,y,x); 





Zero 


(40360) y=(%k2x+%k1)%e É 


Exemplo 2.69 Considere a equação 
v"+2y'+27=0. 


(a) Ache a solução geral; 
(b) Esboce o gráfico da solução que satisfaz as condições iniciais v(0) = 0e v(0)=1. 


Solução: 


Maxima 
ego: diff(y,x,2)+2x"diff(y,x)+2+xy =0; 
2 


d d 
%361) = v+2[>v|]+2y= 
(40861) ay [5x)* nO 
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sol5:ode?2(eqD,y,x); 

(%0362) y=%e *(Y%kl sin(x) + Yk2cos(x)) 
so IP: ie(sols,x=0,y=0, diriiy,s)=1); 

(%0363) y=%e*sin(x) 
yvix):=" Frhs(solPJ: 

(%0364) y(x):=rhs(solP) 
wepilotediyo le So Sly 2] lylabero sos 
[gnuplot preamble,"set size ratio 1; set zeroaxis;"]); 


ho 


15 


y6x9) 


(%0365) 


Exemplo 2.70 Ache a solução geral da equação 


v"+3y'+2y =x. 


Maxima 
eq6: diff(y,x,2)+3x"diff(y,x)+2xy =x; 
E d 
O) ao Noto pi 
RA [5x] V=x 
sol6:0de?(eg6,y,x); 


(%0366) 


DN 
(%0367) y=%k1l%e* + %k2 Poen2R ade 
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Exemplo 2.71 Ache a solução da equação 
v"-2y'+27=0 


que satisfaz as condições de fronteira v(0) = 0 e v(r/2) = 1. Esboce o gráfico da 
solução. 


Solução: 


Maxima 


eq/:"diff(y,x,2)-2x"diff(y,x)+2+xy =0; 
2 d 
%0368) —v-2|] — Dq) 
(60368) 57 (5x) y 


sol7:ode2(eg7,y,x); 


(%0369) y=%e!(Y%k1 sin(x) + %k2 cos(x)) 
SoJPO2: bec(sol7;x=0,94=0 0 Gpii2v=l)s 
r 
(%0370) y= %e o sin(x) 
Ylx):= rhs(50]P02): 
(%0371) y(x):=rhs(solP02) 
wxplotzd(y(v) x, 02] [y5 0,2], Ivy labelo vor 1; 
[gnuplot preamble,"set size ratio 1; set zeroaxis;"]); 


ho 





(%0372) 
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Exemplo 2.72 Encontre a solução y(x) do seguinte problema de valor inicial 


(x-Dy+(y-1)xºy'=0; 
v(3)=2. 


Solução; 


Maxima 


eq8: (x-1)xy23+(y-1)+x23x" diff(y,x)=0; 


d 
(40373) x (p-1)oy+(x-1)p?=0 
sol8:ode?(eg8,y,x); 


Zu = 2x-—1 

Areada ar 

IcllsolB,x=2,)="9); 
27-10 x2+72x-36 


%0375 - 
CD 72%? 


(%0374) 











Exemplo 2.73 Encontre a solução v(x) do seguinte problema de fronteira 


Solução: 


Maxima 
eq dafEly,x,2) E ve diffiy;x)ca = 0; 
d? oo 
%0976) — —v| =0 
PROC av es(ão) 
sol:ode2(eg,y,x); 


3 
6%k1 
ara e oo 
Be2(sol x=0,v=1,x=19453); 
3 
a 
o o e A 


6 2 
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Exemplo 2.74 Calcule y'(x) para v(x) dada implicitamente por vt +x2y? + x! =3. 


Solução: Como a solução y(x) é dada implicitamente, derivaremos a equação usando 
o comando diff e depois resolveremos a equação resultante usando o comando solve. 


Maxima 
(%1454)F:yM+x02xy-2+x04-3=0; 


(%0379) y px sp +x'-3=0 
(%i455)depends(y,x); 


(%0380) [y(x)] 
(%i456) diff(F,x); 


(%0381) DR Do e rr = () 
x 

(%i457)solve(%060,diff(y,x)); 

pe pd 


d 
ROO q 


2.5.1 Campo de direções 


Muita informação importante sobre a equação diferencial de primeira ordem pode 
ser obtida por uma simples análise geométrica da função f(x,y). 

Observe que, em cada ponto do plano (x,y), a inclinação da tangente da solução 
v(x) é dada por f(x,y). O campo de direções é o desenho no plano, de vetores 
tangentes com a inclinação descrita por f(x,y). 


Para desenharmos campos de direção com o Maxima usa-se o pacote plotdf(f(x,y)). 
Exemplo 2.75 Desenhe o campo de direções da equação y' = 1 +y+y?. 


Solução: 


Maxima 
(%i458)1oad(plotdf)$ 
(%i459)plotdf(1 + y + y>2); 
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(%0383) 


= 
y 


Exemplo 2.76 Desenhe o campo de direções da equação y' 


Solução: 


Maxima 


plotdf)$ 
plotdf(-x/y); 


load( 





(%0384) 
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2.5.2 Sistema de equações diferenciais 


Para resolvermos um sistema de equações diferenciais de primeira ordem do tipo 


vi = f(x,y); 
vo = f(x,y); 


Ym = Ím(X 9). 


usamos o seguinte comando: 








desolve([eqnl,eqn2,..., eqnm), [yi(x), y>(x),..., Ym(x)]) 




















A forma de declaração das funções das equações diferenciais para o comando 
desolve devem ser explícita. 


Exemplo 2.77 Encontre a solução do seguinte sistema de equações 


u'=v-zuv=z-u;z'=u-v 
u(0) = 1;v(0) = —1;2(0) = 1. 





(%0389) ala) =u(x)—v(x) 


atvadde (su = ds 
(%0390) 1 


sist: [egni,eqn2,eqn3]; 


d d d 
(%0391) Lg) =v(x)-z(x), Tx (8) =z(x)—u(x), Er) =u(x)—v(x)] 
solSt:desolve(sist,[u(x),v(x),z(x)]); 





(40392) [u(x) = ont E 
= sin(V3x) 2cos(V3x) 1 
a(o) TE + ; E a 


Exemplo 2.78 Encontre a solução do seguinte sistema de equações 


Maxima 
egqni: drfilu(x),x,2)J=v(x); 
2 


(%0393) ua) ESUMN) 


atvalue(u(x),x=0,1); 
(%0394) 1 
atvalue( ditf(u(x),x),x=0,2); 


(%395) 2 
egn2:“diff(v(x),x,2)=u(x); 


2 


(%0396) pie) = (x) 


eo 
atvalue(v(x),x=0,2); 
(%0397) 2 
atvalue('diff(v(x),x),x=0,2); 
(%0398) 2 
Sisto legni eqn2]; 
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2.6 Integração Numérica 


A Integração numérica é um conjunto de técnicas numéricas para o cálculo apro- 
ximado da integral definida da função f. 

Estas técnicas são importantes nos casos, por exemplo, que a função f tem uma 
expressão analítica muito complicada ou é dada por uma tabela. Em geral, uma 
regra de integração numérica consiste em aproximar a integral definida por uma 
soma finita, supondo-se que f é uma função “bem comportada”. 

Para isto, considere uma partição P:a=xg<x/ <..<x;<...<xy = b do intervalo 


de integração [a,b] e a seguinte aproximação: 


b n 
J flcdx=5 wi f(x). (2.2) 
u 1=0 


Os pontos x; são chamados pontos de integração e w; são os pesos da fórmula de 
integração. Cada escolha de x; e w; define uma regra de integração. 

Nesta seção, descreveremos algumas regras clássicas de integração numérica e 
usaremos o Maxima nas implementações. Para maiores detalhes consulte [3]. 


Trapézio repetida 


Considere uma partição P do intervalo [a,b]. Aplicando a regra dos trapézios em 


cada subintervalo [x; 1,x;] obtemos 


fre a=5 (o figdx = 5f(x0) + 2/00) 42500) +... 2f(1) + FO) 


E o erro é dado por 
h2 
Err=5(["(8) ++ 1"(By)) 
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com B; € (x; 1,x;). Usando o teorema do valor médio para integrais tem-se 


h2 
Err =—(b-— 5107), veta, b). 


Exemplo 2.79 Usando a regra dos Trapézios repetida calcule um valor aproximado 
1 


da integral J e dx. 
0 


Solução: 


Maxima 
nodos a: 0 bo o he(b=a) im 
xx:makelist(a+ixh,i,0,n); 


mas os 74 
(Rod O dn [os e as ass 1] 


f(x) :=exp(x-2)$ 
TR:n/2x(sum(2xf(xx[1]),i,2,n)+f(xx[1])+f(xx[n+1])),numer; 


(%0402) 1.467174692738799 
Estimativa do erro: 
Maxima 


define(dTRf(x),abs(diff(f(x),x,2)))$ 
plot2d(dTR£(x),[x,0,1]); 





9 
a cad 
o O 
T 


do cad ad 
= 
T 


po E Ooo om 
——— + ——+—+— 


AXO" XADM 











(%0403) 
ErroTR: (abs(b-a)*xh22x17)/12, numer; 


(%0404) 0.014166666666667 = |ErroTR]| < 0.014166666666667 


13o 


Simpson repetida 


Considere uma partição P do intervalo [a,b]. Para aplicarmos a regra de 
Simpson devemos subdivir o intervalo em um número par de subintervalos, pois 


cada parábola requer três pontos de interpolação. Assim, 


frias = MN Fixidx 


i=1  26-1) 


IR 


(POr) + AF 0a) +2F 00) + 420002) + 4f (1) + FO) 


IR 


h n/2 (n-2)/2 
Brisa rec? Ê. fear) 
i=1 i=1 


E o erro é dado por 
nó e(4) 


Exemplo 2.80 Usando a regra de Simpson repetida calcule um valor aproximado 


1 
da integral H! e“ dx. 
0 


Solução: 





Estimativa do erro: 
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Fórmulas gaussianas 


Fórmulas gaussiana são regras de integração cujos pontos não são igualmente 
(xo,X1,...,Xn) não são necessariamente igualmente espaçados. Serão consideradas 
aqueles exatas para polinômios de grau menor ou iguala 2n+1. 

As fómulas gaussinas mais conhecidas são as chamadas Gauss-Legendre por esta- 


rem associadas aos polinômios de Legendre 


1 
2 
Polx)=1, qnlx)=2, qpalx)=2"—, 
Para obtermos a fórmulas gaussiana de dois pontos vamos considerar as raízes 
1 R 1 1 
de p>(x) = x? — 3º OU seja, xo = 1/7 e X1 = Vi Para obtermos os pesos wq e w 


devemos garantir que a fórmula 
1 
| ftddxewo ft) + ft) 
| 
seja exata para polinômios de grau < 1, ou seja, f(x)=1 e f(x)=x. Assim, 
1 
J ldx=wo+wW1 > w+wW=2 


1 


1 
Wo w 
| vas > wo w/=0. 
- 3 


LS 


Logo, wo=wj=1 ea fórmula Gauss-Legendre para dois pontos é dada por 


Pottvas=s [all] 


que é exata para polinômios de grau < 3. 
Analogamente, fórmula Gauss-Legendre para três pontos é dada por 


ficas (5 /(-N06)+8f(0)+5f(V06)) 


que é exata para polinômios de grau < 5. 
Para usarmos as fórmula gaussinas no intervalo [a, b] fazemos a seguinte mudança 
de variável: 


ie alto =dpsiis a)) 


Assim, 


b-a b+a b-a b-a b+a 
fo f(x ne [s( HE Jax= 2523 E nf ( Ens 3 2) 
Exemplo 2.81 Usando a regra gaussiana calcule um valor aproximado da integral 


l 2 
J E” ar 
0 


Solução: 
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Programas do Maxima para integração numérica 


Podemos também usar os programas do Maxima para calcularmos numericamente 


a integral de uma função. Por exemplo, podemos usar o seguinte comando: 








quad qags(f(x),x,a,b, opções) - calcula numericamente a integral de 
f(x), usando regras do tipo de quadra- 


tura Gaussiana 




















O comando quad qags retorna uma lista de quatro elementos: 
* uma aproximação para o valor da integral; 
* oerro absoluto estimado da aproximação; 
* o número de avaliações do integrando; 
* um código de erro. 


O código de erro são os seguintes: 


O - nenhum problema foi encontrado; 


1 - foram utilizados muitos subintervalos; 


2 - foi detectato um erro de arredondamento excessivo; 


3 - o integrando se comportar muito mal; 


6 - a entrada não é válida. 

As “opções” do comando quad qgags são as seguintes: 
ordem: é um inteiro m que define a ordem da regra de integração; varia em (1,2,3,4,5,6). 
erro: erro relativo desejado; padrão entre -1 à -8; 


particao: número máximo de subintervalos;padrão 200. 
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Para ilustrar o comando quad qags, tentaremos calcular a integral 


1 
J arctg(x)? dx 
0 


com o comando integrate: 





A resposta do aplicativo (retorno da integral definida) indica que o comando in- 
tegrate não é capaz de calcular uma primitiva da função f(x) = arctg?(x). Vamos 
obter um valor aproximado usando o comando quad qgags. 





Observação 2.2 Os comandos quad qapi e romberg também são programas para cal- 


cular numericamente a integral de uma função. 





3.1 Elementos da Geometria Fractal 


Benoit Mandelbrot, matemático polonês, criou este termo em 1975 para designar 
conjuntos com características muito especiais. O termo vem do latim, do adjetivo 
fractus, derivado do verbo frangere que significa quebrar, fracionar. 

Muitos matemáticos referem-se aos fractais como conjuntos de dimensão não in- 
teira, conjuntos de medida Hausdorff, conjuntos com estrutura fina ou simplesmente 
conjuntos irregulares. 

Algumas propriedades caracterizam as figuras que chamaremos de fractais. Sa- 
bemos que uma figura é um conjunto de pontos. O círculo é um conjunto de pontos 
do plano, a esfera é um conjunto de pontos no espaço tridimensional. Estas figuras 
são estudadas na geometria clássica. Porém, os fractais são figuras com proprie- 
dades e características bem peculiares, que os diferenciam das figuras geométricas 
clássicas. Tais figuras são construídas, em geral, por processos iterativos e, seguindo 
Kennet Falconer [5], classificaremos uma figura como um fractal quando apresentar 
todas ou a maiorias das seguintes propriedades: 


Lei de formação simples e complexidade 


O processo de construção de um fractal é geralmente iterativo ou recorrente, po- 
rém relativamente simples. Porém, a simplicidade do processo não é suficiente para 
descre- ver o fractal como um lugar geométrico de pontos que satisfazem uma pro- 
priedade simples ou que são soluções de equações simples. 

Os fractais podem ser construídos por equações simples, mas não podem ser des- 


critos analiticamente por uma equação, como, por exemplo, a reta e a circunferência. 


Estrutura fina em qualquer escala 


A estrutura fina é a propriedade da “complexiadade em escalas”, ou seja, se su- 


cessivas ampliações da figura forem feitas, mais detalhes sobre a mesma serão des- 
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cobertos. 

As figuras geométricas clássicas não possuem esta estrutura. Por exemplo, ao exa- 
minarmos o gráfico da função f(x) = senx no intervalo [0,27], observaremos uma 
série de detalhes tais como pontos de máximo e mínimo, raízes, pontos de inflexão, 
curvatura, etc. Agora, se restringirmos o estudo ao intervalo [0,1], perceberemos 
somente uma ligeira curvatuta e nada mais. Se reduzirmos ainda mais o diâmetro 
do intervalo, digamos a apenas alguns milésimos de radiano, observaremos que o 
gráfico (da senóide) se assemelha ao gráfico de uma reta, e a reta é uma “figura sem 
detalhes”. Portanto, à medida em que a região de investigação é reduzida, o deta- 
lhamento tende a desaparecer e então podemos concluir que a senóide não tem a 


propriedade de estrutura fina. 


Dimensão topológica e dimensão fractal 


As figuras geométricas clássicas têm dimensão fractal igual à dimensão topoló- 
gica, mas os fractais têm, em geral, dimensão fractal estritamente menor que a sua 


dimensão topológica. 


Autossimilaridade ou auto-afinidade 


A autossemelhança ou autossimilaridade é uma ideia antiga e uma propriedade 
geométrica simples. É a “simetria através das escalas”, ou seja, um objeto possui 
autossemelhança se apresenta a mesma forma em qualquer escala em que seja ob- 
servado. Naturalmente, nem todos os objetos geométricos têm esta propriedade. 
Por exemplo, um círculo numa escala muito grande não é nada mais do que uma 
reta. Por outro lado, um quadrado é um conjunto autossimilar do plano, pois pode 
ser formado por quatro cópias deles mesmo reduzidas por um fator r. 

Matematicamente, a propriedade de autossimilaridade é descrita por transforma- 
ções cujas imagens dos objetos sejam “cópias reduzidas” do mesmo e “reconsti- 
tuam” integralmente o mesmo objeto. Tais transformações são chamadas transfor- 
mações similares. 

A auto-afinidade é um tipo de autossimilaridade mais geral. Uma figura é auto- 
afim se pode ser decomposta em partes menores que são não exatamente iguais, 
mas apenas modificadas por transformações afins tais como: contrações, dilatações, 
translações, rotações, reflexões ou combinações dessas transformações. Por exem- 
plo, se uma figura é reduzida uniformemente em todas as direções formará uma 


réplica geometricamente semelhante à figura original. Um círculo que é reduzido 
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deste modo formará um círculo menor. Porém, se a contração se efetuar uniforme- 
mente em uma só direção, resultará em uma elipse. Um retângulo só será similar a 
outro retângulo que possua mesma razão entre o comprimento da base e sua altura, 


enquanto será afim a qualquer paralelogramo. 


3.2 Sistema de funções iteradas 


Nesta seção descreveremos as principais propriedades matemáticas dos fractais. 
Para mais detalhes consulte [4) [5]. 


Definição 3.1 Seja r > O um número real, dizemos que a transformação f: X > X é 


uma “transformação similar” se 








IH-fl="|-7), vzpeR? 
Ser<1, f é chamada uma “contração”. 


Exemplo 3.1 Considere um triângulo equilátero A de lado igual a 1 e a transfor- 
mação que reduz A de tal modo que cada lado seja a metade do anterior. Então, 
f:A->Rº? é dada por 





fosp)=5(,9) 


Observe que, se aplicarmos a transformação f infinitamente temos que os triân- 
gulos serão reduzidos de modo que seus lados serão cada vez menores, ou seja, a 


sequência de triângulos tenderá a um único ponto. 


Figura 3.1: Redução 
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Definição 3.2 Dizemos que dois conjuntos são autossimliares se um é a imagem do outro 
por transformações similares. O conjunto E é autosimilar se é formado pela união de m 


imagens similares de si mesmo, ou seja, 
E=[JfB). (3.1) 


Um conjunto autosimilar definido por (3.1) é chamado atrator ou conjunto invariante 
do sistema de funções iteradas-IFS (f, f>o,..., fu) com fr funções similares com raios 


sro csTm- 


A teoria das IFS (veja, [5]) garante que se o espaço (X,||:||) é completo e [f, fo,..., fm) 
são contrações então existe um único compacto não vazio dado por (3.1). Deste 
modo, temos um processo iterativo de construção do atrator E: dado qualquer com- 
pacto não vazio inicial Kg, para Yn, considere 


Ku = film). 
k=1 


Então, a sequência (K,] converge para E (num certo sentido). Os seguintes resul- 


tados garantem estas afirmações: 





Teorema 3.1 Considere X um espaço métrico completo e K(X) a família de subconjuntos 
compactos de X. Se fi: X — X são contrações com fatores de contração r; e (0,1), 
i=1,2,...,m. Então a aplicação F : K(X) > K(X) definida por 








F(K) =| HE) 
k=1 


tem um único atrator K = F(K). 


O atrator dado no Teorema[3.1|pode ser obtido pelo método das sucessivas aproxi- 
mações. De fato. 





Corolário 3.1 Dado qualquer Kg e K(X). Se 
Kas = F(Km) =) f(Kn) (3.2) 


para n > 0, então a sequência (K,) converge, na métrica Hausdorff, para o atrator do IFS 


Ui )orera fm 
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O Corolário[3.1|justifica o fato que chamarmos o ponto fixo K de F de atrator. 


Agora, também podemos justificar o uso do termo “sistema de funções iteradas”. 





A iteração tratada aqui é, por exemplo, aquela que dado qualquer x, € R” aplicamos 


repetidamente as funções f. em qualquer ordem. Mais especificamente, sejam xy € 





R” e (x,) uma sequência definida por 


Xy = a (44=1) 


paran>1. Então, 

(i) todo ponto de acumulação da sequência (x,) pertence ao atrator K de F; 

(ii) todo ponto do atrator K é um ponto de acumulação de cada sequência (x,) 
para alguma escolha k e (1,2,...,m]; 

(iii) existe um ponto xy e uma escolha de k € (1,2,...,m) tais que K é igual ao 
conjunto de todos os pontos de acumulação de (x,,). 

Uma versão mais sofisticada do resultado (iii) diz que uma escolha “aleatória” 
de k e (1,2,...,m) (com probabilidade 1) terá como fecho o conjunto K (para mais 
detalhes consulte “jogo do caos” em [4)[5]). 

Podemos construir um objeto autossimilar reduzindo (ou ampliando), rodando 
e/ou transladando este objeto. Por exemplo, considere um quadrado O unitário, se 
queremos reduzir este quadrado a metade (fator de redução 1/2) devemos “trans- 
formar” todos os pontos deste quadrado em novos pontos cujas coordenadas serão 
a metade das coordenadas dos pontos originais (veja figura [2.2). Se queremos trans- 
ladar O duas unidades a direita da sua posição inicial devemos “transformar” todos 
os pontos de O em novos pontos cujas coordenadas serão as coordenadas originais 
mais duas unidades. No que segue, descreveremos estas transformações. 

Por exemplo, se P = (x,y) é um ponto de Q, para reduzirmos O a metade transfor- 


maremos P = (x,y) no ponto P, = (x1,y1) tal que x, = 5 evi= - Podemos represen- 


tar esta transformação do seguinte modo: f(x,y) = (5. =) 
O mesmo vale para a translação. Transformaremos o ponto P = (x,y) no ponto P, = 


(x,YV1) tal que x, =x+2ey =. Esta transformação tem a seguinte representação: 
fly) =(x+2,9). 
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tJ| ms 


Figura 3.2: 


Agora, se quisermos rodar O um ângulo 0º no sentido anti-horário, cada ponto 
P = (x,y) de Q deve ser transformado no ponto P, = (x1,y1) tal que x, = x cos0 — 
vsen0 ev; = xsenO + y cos0 (veja figura [3.3). Esta transformação tem a seguinte 
representação: f3(x,y) = (x cos0 —ysen0,xsenO + y cos 0). 


Figura 3.9; 


Observe que as transformações acima podem ser representadas na seguinte forma 
matricial: 


dj 12 x bi 
flwv)= + 
dj 492 v b> 
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De fato, 
é 0 x 
filxy)= É 1 Ê 
1 0 x 2 
folx,v)= + , 


cos0 -senQ x 
falx,v)= 


senô cos0 V 


3.3 Fractais com Maxima 


Nesta seção descreveremos alguns fractais bem conhecidos, os quais implementa- 
remos com o aplicativo Maxima. Estes fractais podem ser implementados usando-se 
os seguintes processos: processo geométrico ou de remoção, processo aleatório ou 


jogo do caos e processo de autossimilares. 


Jogo do caos Escolhe-se os vértices do conjunto da etapa inicial da construção do 
fracatl, re- presentado por ((p1,91),...,(Pm Qm)). Desenha-se o ponto inicial (x9,V0) € 
escolhe-se aleatoriamente um dos pontos de [(p1,91),...,(Pm»Qm)). Então, desenha-se 
um novo ponto que estará no segmento entre o último ponto desenhado e o ponto 
que foi seleccionar aleatoriamente, a uma distância do ponto selecionado. Esta dis- 
tância será 8 vezes o comprimento do segmento, com o fator de redução do fractal. 
Repete-se novamente o processo, um número grande de iterações. 

O pacote dynamics do aplicativo Maxima tem o comando chaosgame que imple- 


menta o jogo do caos. 


Fractais por autossimilaridade Podemos construir fractais usando as transfor- 


mações autossimilares do fractal. Estas transformações são do tipo 


k k k 
41 12 x bi 
frlaop) = k k k 
494 92 y db; 


para k=1,2,3, yin: 
Assim, cada fractal é definido por m matrizes Ap-[ai;) e m vectores by=[b]. E 
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neste caso, o jogo do caos é como segue: escolhido um ponto inicial (x9,v09) aleatori- 
amente, aplica-se, também aleatoriamente, uma transformações autossimilares em 
(x9, Vo) para gerar o ponto (x,,71) e repete-se o processo. Deste modo, na n-ésima ite- 
ração teremos gerado (x,,Y,) pontos do fractal ou muito próximos dos fractal. Como 
o fractal é o conjunto atrator destas transformações, independente do ponto inicial e 
para um número muito grande de iterações, o resultado final será uma figura muito 
próxima ao fractal. 

Para controlar melhor a densidade de pontos em diferentes partes do fractal, con- 
vém que algumas das transformações autossimilares sejam escolhidas com maior 
frequência. Para isso usa-se diferentes probabilidades para a escolha de cada trans- 
formação. 

O comando ifs, no pacote dynamics do aplicativo Maxima, implementa a constru- 
ação de fractais por autossimilaridade. Os parâmetros de entrada do programa são 
uma lista de probabi- lidades, uma lista de matrizes Ay, uma lista de vetores by, um 
ponto inicial (x9,0) e o número de iterações que deverão ser realizadas. 

O pacote fractals também gera alguns fractais por autossimilaridade como 


exemplificaremos nesta seção. 


Programando fractais 
Podemos usar os recursos de programação do aplicativo Maxima e implementar 
fractais, sem necessidade de usarmos pacotes. Apresentaremos alguns exemplos na 


Seção [3.4] 
3.3.1 O conjunto de Cantor ternário 


Construção por remoção 


O conjunto triádico de Cantor é o conjunto formado pela remoção de uma sequên- 
cia de intervalos abertos do intervalo [0,1]. Primeiro, dividimos intervalo [0,1] em 


três partes iguais e removemos a parte do meio, isto é, o intervalo (1/3,2/3) para 


o =[p3Jofpa) 


Novamente removemos o terço médio de cada intervalo fechado de C;, para obter 


Gs=[0.5]u5-5]9[5r5]9 [571] 


obter o seguinte conjunto: 


o seguinte conjunto: 
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Note que C; é a união de 2? = 4 intervalos fechados, cada um do tipo [j/32,(j + 
1)/3?] e que C, 2 C,. 

Novamente removemos o terço médio de cada um destes intervalos para obter 
o conjunto C3, que é a união de 2º = 8 intervalos fechados tal que Cj D C, D Cs. 
Continuando este processo n-estágios obteremos o conjunto C,, que é a união de 2* 
intervalos do tipo [j/3”,(j+1)/3"]talqueC,20,D...DC,D... 

O conjunto triádico de Cantor é a interseção dos conjuntos C,, n EN, ou seja, 


Ee 


neN 


A figura abaixo ilustra o processo de construção do conjunto de Cantor. 


Construção por autossimilaridade 





Vamos aplicar o Teorema Para isto, considere X = R com a métrica usual e 
. . 2 
F: K(M) > K(M) dada por F(K) = A(K)U b(K) com f(x) = 5 e f(x) = 5* 3º 

Ko = [0, 11. 
x 1 xo 2 


Observe que, fi(x) = E efetua uma redução de [0,1] com fator 3 folx) = 3 + 3 


1 
efetua uma redução de [0,1] com fator — e depois uma translaçao de 3 para a direita. 
Fazendo C,,1 = A(C,)U f(C,) temos que 


1 2 1 ad » A 8 
C4=|0,-=|U|=,1|,C,=|0, =]Ul=,=|Ul=,=|Ul=,1],... 
= [osjulga e =[2-5]v[5:5)“[5r5]vl! | 
Observe que f, e f, são transformações similaridades e pelo Teorema F tem 
um único ponto fixo. Mostraremos que este ponto fixo é o conjunto ternário de 











Cantor. 
Para isto, observe que, por indução, obtemos C,., = A(C,)U fo(Cy) para n > 0. 
Primeiro, mostraremos que C € A(C)JU f(C). 
Seja x € C. Então, x € Co e logo, x E Cy, ou seja, x E [o 5| ou x € [5 1) Vamos 
neN 


2 E as E 
considerar x E [51] (o outro caso é análogo). Assim, para todo n sabemos que 
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x EC = h (Cn) U f(Ca). Mas 


1 


filCa E f((0,1])= b 5) 


Logo, x e f>»(C,), ou seja, 3x — 2 € C, para todo n, então 3x — 2 E EG, =C,ou seja, 
neN 
x € fi(C). No outro caso, temos que x e f>(C). Portanto, x e A(C)U f(C). 


Agora, mostraremos que CD A(C)U H(C). Seja x e A(C)U H(C). Então, xe f(C) ou 

x e fb(C). Considere x e f(C) (o outro caso é análogo) então 3x — 2 € C. Assim, para 

qualquer n tem-se 3x -2€C,, ouseja, xe H(C,) CC,. Logo, x E (Cn = E 
neN 


Exemplo 3.2 Desenhe o Conjunto de Cantor usando o jogo do caos 


Solução: 





O aplicativo salvou um arquivo com o nome “Cantor.eps”, no formato EPS, na 


pasta de trabalho. 





contraction factor: 0.33333333333333 





0 01 0.2 03 04 0.s 0.6 07 08 09 1 
The chaos game wrth 2 pomts 


Observação 3.1 Note que existem pontos fora do conjunto. Isto acontece por que o jogo 
do caos é probalilistico. 
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3.3.2 O Triângulo de Sierpiúski 


Construção por remoção 


O processo de construção do Triângulo de Sierpiúski é similar ao da construção 
do conjunto triádico de Cantor: a partir de um objeto inicial, retiramos uma parte 
central sua, e ainda obtendo cópias similares ao objeto inicial. Procedemos então de 
forma similar com cada cópia e assim indefinidamente. Mais precisamente, inicia- 
se com uma região triangular equilátera So, de lado 1 e vértices v; = (0,0), v> = (1,0) 
e vs = (1/2, 3/2), subdividindo-a em quatro regiões triangulares menores, usando 
para isto os segmentos de reta que unem os pontos médios de cada lado. Remove-se 
então o interior da região triangular central (a região central, sem sua fronteira). 
Note que assim é gerado um conjunto S, formado por três regiões triangulares me- 
nores e congruentes, cujos lados medem a metade da medida dos lados do triân- 
gulo original. Repete-se então esse mesmo procedimento em cada uma das três 
regiões triangulares e assim sucessivamente, gerando-se uma sequência de conjun- 
tos SyD S/D S; D-:: com cada S, formado por 3” regiões triangulares, similares à 
região inicial, cujos lados medem 2” da medida do lado do triângulo inicial. 


O Triângulo de Sierpiúski S é o conjunto limite desta sequência de pontos, isto é, 


S=[)Sw 


neN 


A figura abaixo ilustra o processo de construção do Triângulo de Sierpiúski: 


A Pt A À E 
AA bd dd, 


Construção por autossimilaridade 





Para construírmos o triângulo de Sierpiúski por similaridade, considere X = R? 








com a métrica usual e F: K(X) — K(X) dada por F(K) = A(K)U f(K)U fi(K) com 


fi(lsy) = (5.5); 
po = (518) 
aca = (E) 
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e Ko uma região triangular equilátera de lado 1 e vértices v; = (0,0), v> = (1,0) e 
vs =(1/2,N3/2). 


1 
Observe que, f(x,y) reduza Ky a metade, f(x,y) reduza Kg a metade e transla 5 
1 1 
para a direita e f3(x,y) reduza Ky a metade, transla z para a direita e > para cima. 


Fazendo S,,1 = f(S,)U H(S,) U f(S,) temos que o conjunto S, é formado por 
três regiões triangulares menores e congruentes, cujos lados medem a metade da 
medida dos lados do triângulo inicial. Deste modo, gera-se a sequência de conjuntos 
So 2 S&4 D &, D --- Do mesmo modo, como procedemos para o Conjunto de Cantor 
ternário, podemos provar que o ponto fixo de F é o triângulo de Sierpiúski. 


Exemplo 3.3 Desenhe o Triângulo de Sierpisúki usando o jogo do caos. 


Solução: 





O aplicativo salvou um arquivo com o nome “triangSierpinski”, no formato EPS, 
na pasta de trabalho. 


contraction factor: 0.5 








o 04 02 03 04 0.5 06 07 08 D9 1 
The chaos came with 3 points 
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Exemplo 3.4 Desenhe o triângulo de Sienpiúski pelo método de IFS, sabendo que 
as IFS deste fractal são dadas por: m=m, =m> = m3; 


0.5 0 0.0 Us 0.0 
m= E Pi = 5 P2= 5 P3 = 
0 05 0.0 0.0 0.5 
Solução: 
Maxima 
load(dynamics)$ 


m: matrix([0.5, 0.0], [0.0, 0.5]); 
pt: matrix([0.0], [0.0]); 


0,5 1070 
(%0415) 

Dis 

0.0 
(%0416) 

0.0 


p2: matrix([0.5], [0.0]); 


0.5 
(%0417) 


0.0 
p3: matrix([0.0], [0.5]); 


0.0 
(%0418) 


0.5 
RES(Ii2 2], Imp [pis po, pal, [105,0:5]), 1000; 
[gnuplot preamble,"pointsize=0.5"],[gnuplot term,ps], 
[gnuplot out file,"'SierpinskiMaxima.eps"]); 
(%0419) SierpinskiMaxima.eps 
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Iterated Function System of 3 transformations 





3.3.3 A curva de Koch 


Construção por remoção 

A construção da curva de Koch inicia-se com um segmento de reta K9, tomemos 
por exemplo o segmento [0,1], dividindo-o em três partes iguais e substituindo-se o 
o terço médio por dois segmentos adjacentes (tais que estes formem um “triângulo” 
equilátero sem a base), formando então uma linha poligonal K, de quatro segmentos 
consecutivos, em uma espécie de “ponta de estrela”. Repete-se então o processo em 
cada um dos segmentos de K4, obtendo-se assim K, e continua-se indefinidamente 


gerando uma sequência de conjuntos (K9,K4,K»,...), cujo limite é a curva de Koch. 
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Construção por autossimilaridade 


A curva de Koch é o atrator do sistema de funções iteradas dado pelas similarida- 





des 
som = (82) 
cy - (Ec 
Eta = (to é VA cat) 
ao = (028) 


Exemplo 3.5 Desenhe a curva de Kock pelo método de IFS, sabendo que as IFS deste 
fractal são dadas por: 


0.3333 0 0.0 
my = , Pi = , 
0 0.333 0.0 
0.2357 —0.2357 0.3333 
mo = , P2 = , 
0.2357 0.2357 0.0 
—0.2357 —0.2357 0.6666 
mg = » Pa- , 
—0.2357 0.2357 0.0 
0.3333 0 0.6666 
ma = » P4= 
0 0.333 0.0 


Solução: 
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3.3.4 Curvas de Peano 


Em 1878, George Cantor mostrou a existência de uma função bijetora definida no 
intervalo [0,1] e assumindo valores no quadrado [0,1] x [0,1] e a Matemática nunca 
mais foi a mesma. Em outras palavras, Cantor mostrou que existem tantos pontos 
num quadradro quanto num intervalo. Pela primeira vez postulou-se o problema 
da “invariância da dimensão”. Imediatamente os matemáticos da época se pergun- 
taram se tal função poderia ser contínua. A pergunta é relevante, pois se existisse 
tal função contínua teríamos que o intervalo [0,1] e o quadrado [0,1] x [0,1] seriam 
similares, o que foge totalmente da nossa intuição! 

Em 1879, E. Netto mostrou que tal função é necessariamente descontínua e os ma- 
temáticos passaram a buscar uma função contínuas e sobrejetora. Como uma função 
contínua definida em qualquer intervalo assumindo valores no plano é chamada de 


“curva”, a pergunta na mente dos matemáticos da época era: 


Existe uma curva que “passa” por todos os pontos de um quadrado (e o seu 
interior) com área positiva? 


Em 1890, G. Peano respondeu esta pergunta construindo a primeira curva que 
passa por todos os pontos de uma região quadrada. Tais curvas são chamadas curvas 


que preenchem o espaço ou curvas de Peano. 
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Apesar de Peano ter apresentado a primeira curva que preenche o espaço, foi o 
matemático francês D. Hilber quem, em 1891, construiu, pela primeira vez, tais 
curvas por um processo geométrico e recursivo. 

Nesta seção construíremos algums curvas que preenchem o espaço, e apesar de 
tais curvas não serem fractais (o quadrado não é um fractal!) o seu processo de 
construção usa autossimiliaridade. 

Curva de Peano Vamos construir a curva de Peano usando transformações simi- 
lares. Considere um quadrado unitário (e o seu interior) e a sua diagonal, a qual 
chamaremos de gerador. A curva de Peano é construída fazendo-se nove cópias do 


gerador com as seguintes transformações similaridade: 











mm = (53 
sam = [ED DA 
som = (o 

fio) = (5-B)e(a1)=(5+5.-2+1) 
doa = (Split 
a = (Efe ri-t 
go = (Ego 

im = (EB lina! 
co = (EBBe(5o 3408 


Curva de Hilbert 


Agora vamos construir a curva de Hilbert. Para isto, vamos tomar como gerador 
um segmento de reta que une os pontos (0,0) e (1,0). A curva de Hilbert é construída 


fazendo-se quatro cópias do gerador com as seguintes transformações similares: 
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tio = (68 

sim = (Sbpal-(õãs: 

pi = (ba 
tico = (En 


Exemplo 3.6 Desenhe cinco estágios da curva de Hilbert com o pacote fractals 


Solução: 


Maxima 
load(fractals)$ 
n: 100008 
wxplot2d([discrete,hilbertmap(1)]) 





0.8 
0.6 
04 
0.2 
> 0 
-0.2 
-04 
-0.6 


-0.8 
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(%0421) 
wxplot2d([discrete,hilbertmap(2)]) 


1 
0.8 
0.6 
04 
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(%0422) 
wxplot2d([discrete,hilbertmap(3)]) 
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x 


(%0423) 
wxplot2d([discrete,hilbertmap(4)]) 
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(%0424) 
wxplot2d([discrete,hilbertmap(5)]) 


1 
















































































a EBESERcA spa tnts unos cs als sie 
06 PRERESI RR RAS RR ta 
oa Rs RS 
o EEssSErcBEssdais Sic css scalsEica 
Da Rash dora Sicgass pssdcaaACiASES 
pq 
os poCnNaA Pas anna nnaA ones 
"4 -08 -06 -04 -02 0 02 04 06 08 1 


XxX 


(%0425) 
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3.3.5 O conjunto de Julia e Mandelbrot 


O conjunto de Julia é a ilustração muito interessante de como processos aparen- 
temente simples geram conjunto altamente complexos. Funções definidas no plano 
complexo C tal como f(z) = x? + c, com c uma constante, produzem fractais de apa- 
rência bastante “éxotica”. 

Por simplicidade, vamos considerar f :C — C um polinômio de grau n > 2 com 
coeficientes complexos, f(z)=a,2Z" +ap 21 +-..+ao. 

Seja f* a composta de f com ela mesma k vezes, ou seja, fo...of,e f*(z)é a 
k-ésima iterada f(f(...(f(z)...)) de z. Os conjuntos de Julia são definidos em termos 


do comportamento das iteradas f*(z) para k grande. 


Definição 3.3 (Conjunto de Julia) Definimos o conjunto de Julia de f como a fron- 
teira do conjunto K(f)= [ze C; f*(z) » 00). Portanto, J(f) = OK(f). 


O complementar do conjunto de Julia é chamado conjunto de Fatou. 


Observação 3.2 O símbolo -» significa que a sequência f*(z) não tende para oo quando 


k > oo. 


Aqui estamos interessados nos conjuntos de Julia de polinômios que são fractais. 
O estudo da sequência f*(z) para vários valores iniciais z é conhecido como dinâmica 


complexa. Para mais detalhes consulte [5]. 


Exemplo 3.7 Seja f(z) = 2? então f*(z) = 2%. Não é difícil verificar que 


lim f*(z)=0 para I|z|<1; 


k—s00 
lim f*(2) =oo para |z|>1. 
Além disso, para |2| = 1 temos que f*(z) pertence ao círculo |z| = 1, para todo k. 
Portanto, K(f) é o disco unitário |z| < 1 e o conjunto de Julia é a sua fronteita, ou 
seja, o círculo unitário |z| = 1. Naturalmente, neste caso especial, o conjunto de Julia 
não é um fractal. E 


Vamos considerar conjuntos de Julia de polinômios do tipo f.(z) = 2? +c, com c 
uma constante complexa. 


Para definirmos o conjunto de Mandelbrot precisamos da noção de conjunto co- 





nexo. Dizemos que um conjunto 4 Cc R” é conexo se não existem conjuntos abertos 
tais que UUV contém Ae ANU e ANV são não vazios e disjuntos. Intuitivamente, 


ser conexo significa ser uma única “peça”. 
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Definição 3.4 (Conjunto de Mandelbrot) Definimos o conjunto de Mandelbrot M como 


o conjunto dos parâmetros c para os quais o conjunto de Julia de f. é conexo. Então, 
M=(ceC;J(f) é conexo.). 


O conjunto de Mandelbrot contém muita informação sobre a estrutura dos con- 
juntos de Julia. Por exemplo, podemos mostrar que c e M se, e somente se, fF(0) » 
oo (para mais detalhes consulte [5] p. 223]). 


O conjunto de Julia com o Maxima 


O seguinte comando do Maxima gera a imagem de um conjunto de Julia para 
c=x+1y comxey reais: 








julia(x,y, opções) 




















Ao final da tarefa, o aplicativo exibe a mensagem File julia.xpm was created e cria 
um arquivo gráfico no diretório de trabalho, no formato gráfico XPM. Esse arquivo 
gráfico pode ser aberto com qualquer programa gráfico. 


O seguinte exemplo ilustra o uso do comando julia: 


Exemplo 3.8 Vamos gerar o conjunto de Julia de f.(z)= 2? +c com c = —0.55+ 0.6i. 


Solução: 
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Os pontos pretos na figura acima são o conjunto de Julia do Exemplo[3.8| A origem 
encontra-se no centro do quadrado. A região apresentada corresponde a valores re- 
ais e imaginários menores que 1,3 em valor absoluto. Os pontos que não pertencem 
ao conjunto de Julia foram representados com uma cor, que corresponde ao número 
de iterações antes de a sequência se afastar da origem mais do que duas unidades 
(se depois de 40 iterações isso não tivesse acontecido, o ponto continua pintado de 
negro). 

Para aumentarmos o número de iterações, usa-se a opção levels no comando julia. 


Por exemplo, 


Exemplo 3.9 Vamos gerar o conjunto de Julia de f.(z) = 22+c com c = -0.75+0.1i 
com 36 iterações. 
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Podemos usar a opcção radius para reduzir o tamanho da região mostrada no 


gráfico: 


Exemplo 3.10 Vamos gerar o conjunto de Julia de f.(z)= 22 +c com c=-0.75+0.1i 
com 48 iterações e raio 1. 








Como 160 iterações (veja figura abaixo), eliminamos ainda mais pontos que não 


pertencem ao conjunto de Julia: 
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Podemos gerar conjuntos de Julia com o pacote fractals do seguinte modo: 





O conjunto de Mandelbrot com o Maxima 


O seguinte comando do aplicativo Maxima gera a imagem de um conjunto de Man- 
delbrot: 








mandelbrot(opções) 
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Ao final da tarefa, o aplicativo exibe a mensagem File mandelbrot.xpm was created 
e cria um arquivo gráfico no diretório de trabalho, no formato gráfico XPM. Esse 
arquivo gráfico pode ser aberto com qualquer programa gráfico. 


O seguinte exemplo ilustra o uso do comando mandelbrot: 


Exemplo 3.11 Vamos gerar um conjunto de Mandelbrot com 30 iterações. 





O conjunto de Mandelbrot reune num único gráfico todas as formas encontradas 


nos conjuntos de Julia. Por exemplo, 


Exemplo 3.12 Vamos gerar um conjunto de Mandelbrot com a ampliação de uma 
pequena região do conjunto, centrada no ponto 0.3 + 0.51. 
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Podemos gerar conjuntos de Mandelbrot com o pacote fractals do seguinte modo: 


Maxima 

(%i528)1oad(fractals)$ 

(%i529)n: 100008 

(%i530)plot3d (mandelbrot set, [x, -2.5, 1], [y, -1.5, 1.5], 
[gnuplot preamble, "set view map"], 
[gnuplot pm3d, true], [grid, 150, 150]); 
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3.3.6 Mais exemplos 
Nesta seção, usaremos o pacote fractals para gerar outros fractais. 


Exemplo 3.13 Desenhe o fractal gerado por cinco pontos, nos vértices e no centro 
de um quadrado, com factor de redução 1/3. 


Solução: 
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Exemplo 3.14 Desenhe o fractal “árvore” dado pelas IFS: 


—0.550 —0.179 —(0.438 
m = » Pi=— , 

—0.179 0.550 0.382 

—(0.246 0.193 —0.379 
mo = » P2= , 


0.275 0.365 0.538 
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0.006  —0.014 —0.283 
—0.147 —0.459 0.490 
Solução: 
Maxima 
load(dynamics)$ 
mi: matrix([-0.550, -0.179], [-0.179, 0.550]); 
pt: matrix([-0.438], [0.382]); 
—0.550 —0.179 
(%0434) 
—0.179 0.550 
—(0.438 
(%0435) 
0.382 
m2: matrix([-0.246, 0.193], [0.275, 0.365]); 
pe: matrix(|-0.979]5 10538]: 
—0.246 0.193 
(%0436) 
(0275 0,365 
—0.379 
(%0437) 
0.538 
m3: matrix([0.006, -0.014], [-0.147, -0.459]); 
p3: matrix([-0.283], [0.490]); 
0.006  —0.014 
(%0438) 
—0.147 —0.459 
—0.283 
(%0439) 


0.490 


-0.55 -05 


-045 04 “035 -03 0.25 02 


Exemplo 3.15 Desenhe o fractal dado pelas IFS: 


Solução: 


0.0 0.577 

, Pp1 E 
—0.577 0.0 
0.0 0.577 

, P2 e 
—0.577 0.0 
0.0 0.577 

E 
—0.577 0.0 
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-0.15 -01 -0.05 


0.0951 


0.5893 


0.4413 


0.7893 


0.0952 


0.9893 





0.0 (577 
(%0441) 

-0.577 0.0 

0.0951 
(%0442) 

(0,5893 


m2: matrix([0.0, 0.577], [-0.577, 0.0]); 
p2: matrix([0.4413], [0.7893]; 


0.0 0.577 
(%0443) 

—0.577 0.0 

0.4413 
(%0444) 

0.7893 


m3: matrix([0.0, 0.577], [-0.577, 0.0]); 
po: matrix(]0:0952], [0,9898]): 


0.0 0.577 
(%0445) 

—0.577 0.0 

0.0952 
(%0446) 

0.9893 


TEs([1,2,8],[mi;m2,m5],[p1,p2,p3],[0.0952,0,9893],30000, 
[gnuplot preamble,"pointsize=0.5"],[gnuplot term,ps], 


[gnuplot out file,"dragao.eps"]); 
(%0447) dragao.eps 
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1 
0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 


0.1 





0 
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 


Iterated Function System of 3 transformations 


3.4 Programando fractais com o Maxima 


Também podemos gerar fractais programando como o aplicativo Maxima. Nesta 


seção apresentaremos um algoritmo para gerar fractais pelo processo de autossimilari- 


dade. O algoritmo é o seguinte: 


1. 


ds 


Inserir as funções iteradas; 


Definir uma lista de pontos iniciais; 


. Inserir o número de etapas (ou iterações) que o Maxima deverá executar; 
. Gerar a lista de pontos que compõe uma determinada etapa do fractal; 


- Ligar os pontos por segmentos de reta. 


Exemplo 3.16 Faça um programa para gerar a Curva de Koch. 


Solução: No apêndice[A.1]apresentaremos a programação passo a passo. 
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£[4](X):=1/3=[X[1],X[2]]+[2/3,0]$ 
Koch(n):=(k:n,L[0]:[[1,0]], 
(for m:1 thru k do 

L[m]:create Jist(f[il](L[m-1][51]),i,1,4,5,1,4-(m-1))), 
wxplot2d([discrete,append([[0,0]],L[k1])1,[x,-0.1,1.1], 
[y,-0.1,0.5], [ylabel, "Curva de Koch"]))$ 














Koch(5) 
05 
04 | 
S 03 
g 02 | | 
o | | 
-0.1 
0 02 04 06 08 1 
X 
(%0448) 


Exemplo 3.17 Faça um programa para gerar cinco iterações da Curva de Peano. 


Solução: 


Maxima 
£[1](X):=1/3=[X[1],X[2]]8$ 
£[2](X):=1/3x[-X[2],X[1]1+[1/3,0]$ 
£[3](X):=1/3=[X[1],X[2]1+[1/3,1/3]8 
£[4](X):=1/3x[X[2],-X[1]]+[2/3,1/3]8$ 
f[5](X):=1/3x[-X[1],X[2]]+[2/3,0]$ 
f[6](X):=1/3x[X[2],-X[1]]+[1/3,0]8$ 
£[7](X):=1/3*[X[1],X[2]]+[1/3,-1/3]$ 
£[8](X):=1/3x[-X[2],X[1]]+[2/3,-1/3]$ 
f[9](X):=1/3=[X[1],X[2]]+[2/3,0]$ 
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Maxima 


(%0449) 


(%0450) 


Peano(n):=(k:n LO: [[1,0]], 
(For mal the k do 

Limjcreate Listltji Ln Glosa, 9,1 1,9-(m-lh)), 
wxplot2d([discrete,append([[0,0]],L[k])1,[x,-0.1,1.1], 


[y,-0.5,0.5], [xlabel,"Curva de Peano"],[ylabel,""]))$ 
Peano(1) 





04 





Curva de Peano 
o 





0 0.2 04 0.6 0.8 1 
iteração: 1 


Peano(2) 











Curva de Peano 
o 





0 0.2 04 0.6 0.8 1 


iteração: 2 


Peano(3) 





Curva de Peano 
o 








0 0.2 04 0.6 0.8 1 


iteração: 3 
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(%0451) 
Peano(4) 





Curva de Peano 
[em | 














0 0.2 04 0.6 0.8 1 


iteração: 4 


(%0452) 
Peano(5) 





Curva de Peano 
o 











0 0.2 04 0.6 0.8 1 


iteração: 5 


(%0453) 
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Exemplo 3.18 Neste exemplo programamos um fractal de nossa autoria, o qual cha- 
mamos “Árvore de Natal.” 


Solução: Vamos gerar o fractal “Árvore de Natal”. 


Maxima 
f[1](X):=2/3x[X[1],X[2]]+[0,1]$ 
f[2](X):=2/3x[float(cos(2x%pi/3))xX[1]-float(sin(2x%pi/3))+X[2], 


float(sin(2+x%pi/3))*X[1]+float(cos(2x%pi/3))xX[2]]+[0,1]$ 
f[3](X):=2/3x[float(cos(4+*%pi/3))xX[1]-float(sin(4x%pi/3))+X[2], 


float(sin(4x%pi/3))*X[1]+float(cos(4x%pi/3))xX[2]]+[0,1]$ 
fractalnovo(n):=(k:n,L[0]:[[0,1]], 


(for m:1 thru k do 
L[mJ:create list(f[i](L[m-1][5]),1,1,3,5,1,30(m-1))), 
wxplot2d([discrete,append([[0,0]],L[k])1],[x,-2,2], 


[ylabel, "árvore de natal"]))$ 
fractalnovo(7); 





árvore de natal 
ER 
on 











(%0454) 


A PROGRAMAÇÃO 


A.l Fractais 


Curva de Koch 


Maxima % Inserindo as IES% 
E[1](X):= 1/3xX$ 
£[2](X):= 1/3x[X[1]*0.5-float(sgrt(3))xX[2]x0.5, 


float(sgrt(3))*0.5+X[1]+0.5+X[2]]+[1/3,0]$ 
FLS] (XxX) = 1/32]X]1]+0.5+Float(sqrt(3))24[2]+0.5, 


-float(sgrt(3))*0.5xX[1]+0.5xX[2]]+[1/2, 


float(sgrt(3))/6]$ 
EJA J(X)=1/8+=[X 1 ,X]2]+[2/3,0]$ 


Y% Pontos iniciais Y% 
L[0]:[[1,0]]8 
Y Iterações Yo 
k = 58 
Y%Lista de pontos% 
for m:1 thru k do L[m]:create list(f[i](L[m-1][j]), 


i,1,4,5,1,4>(m-1))S 
Y%Ligando os pontos por segmentos de reta: construção do fractal% 


plot2d( [discrete,append([[0,0]],L[k])],[x,-0.1,1.1], 
1y,=0.1,0.5], [ylabel; “Curva de Koch" |); 
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